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Introduction 

The  decline  in  the  water  tables  of  the  Great  Plains, 
coupled  with  the  rapid  increase  of  irrigated  area  in  this 
region,  has  made  conservative  use  of  water  on  agricultural 
lands  an  imperative.  High  investment  costs  for  water  and 
energy  have  also  made  modern  irrigation  scheduling  techniques 
attractive  to  today's  farmer.  The  soil-water  reservoir  of  the 
Great  Plains  is  primarily  depleted  by  evapotr anspiration  (Et) 
(Kanemasu  et  al.,  1976),  defined  as  the  evaporation  of  water 
from  soil  and  plant  surfaces  combined  with  the  transpiration  of 
water  through  plant  tissues.  Accurate  predictions  of  crop  Et 
throughout  the  growing  season  can  provide  the  irrigator  with 
the  timing  and  amount  of  water  which  must  be  applied  to  prevent 
yield-reducing   stress. 

Basal  crop  coefficient  curves  for  corn  ( Zea  mays  L.)  and 
soybean  ( Glycine  max  (L.)  Merr.  )  were  developed  from  water  use 
and  climatological  data  collected  at  sites  at  Manhattan  and 
Tribune,  Kansas  from  1974  through  1982.  A  crop  coefficient  is 
defined  as  the  ratio  of  crop  evapotranspiration  (Etc)  to  some 
reference  crop  Et  (Etr).  In  this  study,  Etr  values  were 
calculated  by  using  equations  requiring  meteorological  input 
data.  Since        rainfall      events      increase      the        evaporation 

component  of  Etc,  and  soil  water  depletion  tends  to  decrease 
the  transpiration  component,  time  periods  with  precipitation 
amounts      less    than   3   mm  were    used   in   the   model    development      and 


the    crop    coefficients   were    corrected    for    available    soil      water. 

A  plot  of  crop  coefficients  for  a  particular  crop  as  a 
function  of  time  is  called  a  crop  coefficient  curve.  The  time 
scale  commonly  used  is  day3  post  emergence.  Such  curves  must 
be  adjusted  for  different  planting  dates,  locations,  years,  and 
cultivars  of  different  growing  season  lengths.  This  is  due  to 
the  fact  that  climate  and  cultivar  differences  result  in 
differences  in  crop  development,  thus  causing  shifts  in  days 
post  emergence-based  crop  coefficient  curves.  Therefore,  one 
goal  of  this  project  was  to  develop  curves  for  corn  and  soybean 
which  could  be  applied  to  a  variety  of  cropping  situations 
without    the    need    to   adjust    the    curves   for    crop   development. 

The  basal  crop  coefficients  described  above  were  plotted 
against  fraction  of  growing  season  calculated  by  using  heat 
unit      methods.  This      is      the      heat      units      accumulated      from 

emergence  to  the  time  in  question,  divided  by  total  heat  units 
accumulated      from      emergence      to      physiological      maturity.  A 

photothermal  unit  equation  was  used  to  calculate  fraction  of 
growing  season  for  soybean.  The  predictive  ability  of  these 
basal  crop  coefficient  curves  was  tested  by  using  data 
collected   at    the   Manhattan    site    in    1983    and    1984. 

In  addition  to  evaluating  the  validity  of  expressing  crop 
coefficient  curves  based  on  fraction  of  growing  season 
(calculated  by  using  heat  units),  other  goals  of  the  project 
were  to  evaluate  various  heat  unit  methods,  develop  curves  of 
maximum  depth  of  depletion  vs.  fraction  of  growing  season  as  a 
means    of    estimating   effective   rooting   depth,     and    to    develop   and 


test    curves   of    relative   leaf    area    index  vs.     fraction   of    growing 
season   for    use    in   net   radiation   equations. 


LITERATURE  REVIEW 

Crop.  Coefficient  Curves 

Crop  coefficients  are  ratios  of  actual  crop 
evapotranspi ra t ion  (Et)  to  the  potential  Et  of  a  well  watered 
reference   crop.    A   common  form  of  the  crop   coefficient   is: 

Kc  =  Etc/Etr 
in  which  Kc  is  the  dimens i onl e s s  crop  coefficient  for  a 
particular  soil  moisture  condition.  Etc  is  crop  Et,  and  Etr  is 
the  daily  reference  Et.  Theoretically,  Etr  represents  the 
evaporative  demand  imposed  by  the  climate  and  Kc  is  a  measure 
of  the  ability  of  the  crop-soil  system  to  meet  the  demand 
(Wright,  1981).  A  plot  of  crop  coefficients  for  a  particular 
crop  as  a  function  of  time  is  called  a  crop  coefficient  curve. 
The  purpose  of  the  curve  is  to  obtain  an  estimate  of  crop  Et  by 
selecting  the  Kc  specific  to  the  time  period  in  question,  and 
then  multiplying  it  by  an  estimate  of  Etr  calculated  from 
meteorological  data.  The  time  base  and  method  of  calculating 
Etr  must  be  compatible  with  the  methods  used  to  develop  the 
particular  curve  being  used.  Such  estimates  of  crop 
evapotranspira t i on  are  essential  for  accurate  irrigation 
schedul ing . 

Denmead  and  Shaw  (19S9)  discussed  a  method  of  estimating 
water  use  by  crops  which  was  originally  suggested  by  Penman 
(1956).    The   method  involves  multiplying  estimated  free  water 


evaporation  by  an  empirical  factor.  This  factor  can  only  be 
applied  to  a  green  crop,  which,  is  actively  growing,  completely 
shading  the  ground,  never  short  of  water,  and  of  uniform 
height.  In  an  attempt  to  determine  the  period  dnring  the  Iowa 
growing  season  when  these  conditions  are  met,  they  plotted 
ratios  of  measured  corn  ( Zea  m  a  y.  s  L.)  Et  to  open  pan 
evaporation  against  time.  After  planting,  their  curve  assumed 
the  form  of  a  sigmoid  growth  curve  as  leaf  area  developed  and 
more  bare  soil  was  shaded.  At  silking,  leaf  expansion  was 
complete  and  the  curve  remained  flat  for  16  days  at  a  ratio  of 
0.81.  They  concluded  that  this  is  the  period  when  Penman's 
method  can  be  applied  to  corn  in  Iowa.  Shaw  (1963)  later  used 
this  curve  in  a  technique  for  predicting  soil  moisture  changes 
snder  corn.  Weekly  averages  of  daily  pan  evaporation  were 
multiplied  by  the  appropriate  ratio  to  obtain  an  estimate  of 
corn  Et.  Doss  et.  al.  (1962)  also  plotted  ratios  of  Et  to  open 
pan  evaporation  as  a  function  of  time,  thus  developing  a  curve 
for  irrigated  corn  in  Alabama.  The  Et/pan  evaporation  ratio 
was  0.38  at  emergence  and  increased  to  a  maximum  of  1.12  during 
early  dough  stage.  It  then  dropped  to  0.95  at  grain  maturity. 
This  contrasts  with  the  Denmead  and  Shaw  curve,  in  which  the 
Et/pan  evaporation  ratio  began  to  decline  after  silking  and 
commencement  of  ear  growth,  and  ultimately  reached  its  original 
low  level  as  physiological  activity  of  the  plant  stopped. 

Jensen  et  al.  (1970)  recommended  a  crop  coefficient  which 
was  not  based  solely  on  stage  of  growth.  This  crop  coefficient 
was   also  based  on  time  since  an  irrigation  or  rainfall  and   on 


the  remaining  available  soil  moisture.  Doss  et  al.  (1962) 
reported  that  the  Et  rate  of  corn  decreased  with  decreasing 
available  soil  water.  Znr  et  al.  (1983)  found  that  linear 
relationships  between  the  transpiration  rate  of  soybean 
[Glycine  max  (L.)  Kerr.]  and  both  radiation  flux  and  vapor 
pressure  deficit  existed  only  in  well  watered  plants  having  low 
stomatal  resistances.  However,  the  transpiration  rate  quickly 
decreased  following  stomatal  closure  and  remained  essentially 
constant  at  high  stomatal  resistances.  Conversely,  Et  rate 
will  increase  after  rain  or  irrigation  due  to  increased 
evaporation  from  the  soil  surface  (Wright,  1981).  Wright 
(1982)  published  a  crop  coefficient  curve  for  snap  beans 
(Phase  plus  vul.«_a_r  is.  L.)  in  which  the  Kc  values  calculated  from 
lysimeter  data  and  computed  reference  Et  were  plotted  against 
time.  This  curve  is  so  erratic  that  it  would  hardly  be  useful; 
the  Kc  values  rise  and  drop  sharply  after  each  irrigation  and 
rainfall.  Wright  remedied  this  situation  by  developing  a 
"generalized  basal  crop  coefficient  curve".  This  was 
accomplished  by  manually  fitting  a  curve  to  the  time 
distribution  of  Kc  representing  conditions  when  the  soil 
surface  was  dry  and  the  availability  of  soil  water  did  not 
limit  plant  growth  or  transpiration.  The  basal  crop 
coefficient,  designated  as  Kcb,  was  usually  set  equal  to  Kc 
after  full  cover.  When  ripening  began,  wet  soil  effects  were 
again  considered.  Wright  found  that  the  maximum  values  of  Kcb, 
corresponding  to  effective  full  cover,  usually  occurred  shortly 


after  the  rows  closed  for  most  of  the  crops.  This  was  when 
LAI  reached  2.5  to  3.0.  The  Kcb  then  usually  declined  with 
time  becanse  of  plant  lodging  and  natural  senescence. 

The  USDA  -  Agricultural  Research  Service  Irrigation 
Scheduling  Program  described  by  Jensen  et  al.  (1971)  contains 
equations  for  estimating  Kc  from  a  basal  crop  coefficient  and 
correction  factors  for  available  soil  water  and  surface 
wetting.   The  equation  given  is: 

Kc  =  KcoKa  +  Ks 
in  which  Kco  is  the  mean  crop  coefficient  based  on  experimental 
data  where  soil  moisture  was  not  limiting  and  normal  irrigation 
stands  were  used  (it  is  therefore  identical  to  Kcb  described  by 
Wright,  1982).  Ka  is  a  coefficient  related  to  a  natural 
logarithmic  function  of  available  water,  and  is  calculated  as 
Ka  ■  ln( AM+1) / lnlOl  in  which  AM  is  the  percentage  of  remaining 
available  soil  water.  Ka  is  therefore  equal  to  1  when  AM  is 
100%  and  goes  to  0  as  AM  goes  to  0.  Ks  corrects  for  the 
increase  in  the  coefficient  when  the  soil  surface  is  wetted  by 
irrigation  or  rainfall.  It  is  calculated  as  Ks  =  (Kl  -  Kci)e~ 
,  in  which  Kl  is  usually  equal  to  0.9,  Kci  is  the  basal  crop 
coefficient  at  the  time  of  irrigation  or  rainfall,  t  is  the 
number  of  days  after  irrigation  or  rainfall,  and  e  represents 
the  combined  effects  of  soil  characteristics  and  evaporative 
demand.  Wright  (1981)  suggested  the  following  formula  for 
correcting  the  Kcb  value  for  surface  evaporation: 

Kc  =  Kcb  +  (1  -  Kcb)(l  -  ( t / t d) lf 2) ( fw) 
in  which  t  is  the  number  of  days  after  rain  or   irrigation,   td 


is  the  usual  number  of  days  for  the  soil  suface  to  dry,   and  fw 
is  the  relative  portion  of  the  soil  surface  originally  wetted. 

Shaw  (1963)  corrected  Denmead  and  Shaw's  (1959)  crop 
coefficients  by  multiplying  them  by  a  stress  factor.  These 
stress  factors  were  obtained  from  curves  of  Etc/pan  evaporation 
ratios  as  a  function  of  percent  available  soil  water.  Three 
separate  curves  were  derived  for  low,  average,  and  high  stress 
conditions.  These  stress  conditions  represented  the 
evaporative  demand  imposed  by  atmospheric  conditions  and  were 
defined  by  various  levels  of  pan  evaporation.  A  separate  set 
of  curves  was  used  after  31  July  since  root  penetration  had 
presumably  stopped  and  atmospheric  stress  would  therefore 
result  in  greater  reductions  in  evapotranspirat ion.  The 
results  of  a  lysimeter  study  of  evaporation  from  a  corn  canopy, 
conducted  by  Ritchie  (1973),  indicated  that  the  amount  of 
available  soil  water  did  not  influence  evaporation  rates  nearly 
as  much  as  Shaw's  curves  imply.  Ritchie  found  instead  that 
evaporation  rates  were  independent  of  soil  water  status  until 
soil  water  was  depleted  beyond  a  critical  threshold.  From  that 
point,  water  was  extracted  at  a  decreasing  rate  before 
evaporation   practically   stopped. 


The  available  soil  water  percentage  used  to  calculate  Ka 
should  apply  to  the  zone  of  active  root  development.  Several 
attempts  to  describe  the  depth  progression  of  corn  and  soybean 
roots  have  been  made,  some  with  the  purpose  of  '.  defining 
available  soil  water  holding  capacity. 


Foth  (1962)  reported  the  results  of  root  dry  weight 
measurements  of  vertical  soil  samples  of  corn  taken  at  various 
growth  stages.  He  found  that  root  growth  below  38  cm  largely 
occurred  during  the  stages  of  rapid  stem  elongation,  tasseling, 
silking,  and  pollination.  By  early  milk  stage,  the  increase  in 
root  weight  between  38  and  91  cm  had  virtually  stopped.  Shaw 
(1963)  presented  a  method  of  estimating  soil  moisture  under 
corn  in  Iowa  in  which  rooting  depth  was  considered  to  be  as 
follows:  to  27  June,  61  cm;  to  4  July,  76  cm;  to  11  July,  91 
cm;  to  18  July,  107  cm;  to  25  July,  122  cm;  to  1  August,  152 
cm.  Shaw,  therefore,  considered  root  zone  development  to  be  a 
linear  function  of  time  with  the  deepest  root  depth  reached  on 
1  August.  Shaw  used  a  plot  of  the  ratios  of  corn  Et  to  open 
pan  evaporation  as  a  source  for  estimating  Et.  This  curve,  as 
reported  in  an  earlier  article  by  Denmead  and  Shaw  (1959), 
indicated  that  the  average  date  of  silking  in  Iowa  is  20  July. 
The  date  of  1  August  would  therefore  be  approximately  12  days 
after  silking.  Stegman  et  al.  (1977)  computed  Et/Etp  ratios 
for  successive  15-cm  increments  of  profile  depth  under  various 
crops  in  an  attempt  to  estimate  approximate  rooting  depth  with 
time.  The  curve  thus  derived  for  corn  appears  to  level  off  at 
70  days  post  emergence.  According  to  their  crop  coefficient 
curve  for  corn,  this  is  approximately  8  days  past  silking. 
Allmaras  et  al.  (1975)  suggested  that  determinations  of  water 
sink  alone  cannot  predict  rooting  depth,  since  in  their  study, 
the   rooting  depth  of  corn  was  almost  always  greater   than   the 
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maximum  depth  of  water  uptake.  However,  since  the  available 
soil  water  correction  discussed  above  applies  to  the  effective 
root  zone,  maximum  depth  of  water  uptake  would  be  an  adequate 
estimate  of  rooting  depth  for  this  purpose. 

Mitchell  and  Russell  (1971)  reported  the  dry  weights  for 
roots  contained  in  soil  monoliths  taken  from  soybean  rows  at 
various  times  during  the  growing  season.  They  divided  root 
development  into  three  general  phases  related  to  top  growth. 
Dnring  phase  2,  which  included  flowering  and  the  beginning  of 
pod  formation,  lateral  roots  penetrated  the  soil  to  a  depth  of 
122  cm  and  root  dry  weight  in  the  92  to  122  cm  layer  was  0.06 
grams/plant.  During  the  third  phase,  which  included  seed  set 
and  maturity,  major  lateral  roots  elongated  rapidly  to  a  depth 
of  122  to  183  cm  and  root  dry  weight  in  the  122  to  183  cm  layer 
was  0.03  grams/plant.  These  results  are  similar  to  those  of 
Winter  and  Pendleton  (1968).  Their  excavation  of  two  soybean 
plants  grown  on  a  fairly  dry  silt  loam  revealed  11%  of  total 
roots  in  the  91  to  122  cm  layer  and  2%  in  the  122  to  152  cm 
layer.  Hayaki  et  al.  (1976)  also  studied  soybean  root 
development  using  field  samples  of  soybean  root  systems.  They 
found  that  the  soybean  roots  in  both  irrigated  and  non- 
irrigated  plots  reached  160  cm  by  stage  R9  (full  size  green 
beans  in  one  of  the  four  uppermost  pods).  However,  there  was  a 
decrease  in  the  rate  of  root  depth  increase  near  stage  R5 . 5 
(beginning  of  pod  development).  The  soybean  root  zone  depth 
vs.  days  post  emergence  curve  ,  derived  by  Stegman  et  al. 
(1977)  as  described  previously  for  corn,   seems  to  level  off  at 
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80  days  post  emergence  which  they  report  as  upper  pod  fill. 
However,  the  slope  of  the  curve  up  to  this  point  is  quite  low, 
and  only  5  to  8  cm  in  soil  depth  are  gained  between  60  and  80 
days  post  emergence.  This  seems  to  agree  with  the  conclusion 
of  Mayaki  et  al.  (1976)  that  the  growth  of  soybean  roots  slows 
down  considerably  after  beginning  pod  development.  Likewise, 
Mitchell  and  Russell  (1971)  found  only  half  the  accumulated  dry 
weight  at  the  next  lower  depth  after  beginning  pod  development 
(during  phase  3) . 

The  pattern  described  above  is  not  evident  in  a  rhizotron 
study  of  7  soybean  cultivars  conducted  by  Kaspar  et  al.  (1978). 
They  found  that  the  cultivars  differed  in  the  rates  at  which 
their  roots  extended  downward  during  the  various  stages  of 
reproduction  and  in  the  pattern  of  depth  increases  displayed  by 
their  root  systems.  The  roots  of  all  seven  cultivars  reached 
the  bottom  of  the  217  cm  deep  compartment,  one  by  the  end  of 
pod  development,  four  by  the  time  of  beginning  bean  fill,  and 
two  reached  217  cm  after  bean  development  was  nearly  completed. 

Several  predictive  simulations  of  soybean  rooting  depth 
have  been  developed.  Stone  et  al.  (1983)  developed  a 
simulation  based  on  equations  requiring  time  in  days  and  soil 
temperature.  These  equations  were  derived  from  data  obtained 
from  a  greenhouse  study  of  the  root  growth  of  four  cultivars 
grown  at  four  temperatures  (Stone  and  Taylor,  1983).  Narda  and 
Curry  (1981)  developed  a  model  of  soybean  root  growth  and  water 
uptake   based   on   root   growth  attributes  and   the   amount   of 
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carbohydrate  available  for  root  growth.  The  calculated  water 
uptake  rates  were  modified  by  accounting  for  decreased  rates  of 
water  nse  due  to  diminishing  soil  water  content,  as  well  as 
aging  and  death  of  roots. 

Various  time  scales  have  been  used  for  crop  coefficient 
curves  and  their  disadvantages  have  been  noted  by  researchers. 
Stegman  et  al.  (1977)  presented  curves  for  six  crops  based  on 
days  post  emergence.  They  pointed  out  that  the  position  of 
these  curves  can  shift  from  year  to  year  due  to  different 
rates  of  crop  development.  They  suggested  that  technicians 
visit  fields  periodically  to  carefully  observe  phenologic 
stages.  If  needed,  curve  adjustments  could  be  made  by 
comparing  growth  stages  in  particular  seasons  with  the  average 
phenology  vs.  days  post  emergence  relationships  shown  on  their 
curves.  In  an  earlier  paper  outlining  the  use  of  crop 
coefficient  curves  for  irrigation  scheduling,  Stegman  and  Valer 
(1972)  said  that  the  time  scale,  days  post  emergence,  was 
chosen  to  purposely  specialize  the  curves  for  east-central 
North  Dakota.  Wright  and  Jensen  (1978)  published  a  crop 
coefficient  curve  for  beans  in  which  Ec  is  a  function  of 
percentage  of  time  from  planting  until  full  cover.  After  full 
cover,  Ec  is  a  function  of  elapsed  days.  This  was  done  because 
large  differences  in  planting  dates  usually  have  little  effect 
on  the  date  full  cover  is  reached.  Wright  and  Jensen  (1978) 
suggested  that  service  groups  use  some  key  growth  stages  to 
adjust  crop  development  to  this  curve. 

Dylla  et  al.   (1980)  recommended  that  crop  coefficients  be 
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developed  from  data  in  climates  similar  to  those  in  which  they 
are  to  be  used.  They  warned  users  that  even  if  this  rule  is 
followed,  year  to  year  climate  aberrations,  crop  v'arietal  and 
maturity  differences,  and  any  other  factors  that  inhibit  crop 
growth  could  cause  shifts  in  the  calculated  Kc  curve.  Wright 
(1982)  commented  that  it  would  be  desirable  to  have  a  means  of 
relating  crop  coefficients  more  directly  to  crop  development 
with  an  index  such  as  accumulated  growing  degree  days.  In 
1976,  Stegman  and  Olson  developed  a  crop  coefficient  curve  for 
pinto  beans  (Phaseolns  vnl uri s  L.)  based  on  accumulated 
growing  degree  units.  They  felt  that  since  this  parameter  is 
more  specific  to  the  temperature  of  a  given  season,  its  use 
would  reduce  the  seasonal  variation  of  the  curve.  Sammis  et 
al.  (1985)  calculated  crop  coefficients  on  a  monthly  time 
period  and  regressed  these  against  monthly  cumulative  modified 
growing  degree  days  (see  next  section).  They  called  this  the 
6  method.  Their  crop  curves  for  corn  were  statistically  the 
same  over  2  years  and  locations.  They  reported  that  corn 
required  1680  growing  degree  days  to  mature  at  both  locations, 
so  varieties  of  different  growing  season  lengths  were  not  used 
in  their  analysis.  They  also  reported  that  some  variability 
between  years  occurred  that  could  not  be  accounted  for,  showing 
the  "inadequacy  of  the  G  method  to  account  for  all  the 
variability  between  years." 


ES.£erenc_e  Ctoj  Evapotranspirati  on 

As   many  as  50  methods  or  variations  have   been   advanced 
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for  estimating  potential  Et  (Hill  et  al.  1983).  Dne  to  the 
ambiguities  in  describing  potential  Et,  the  term  "reference 
crop  Et  "  (Etr)  is  now  used  by  most  researchers,  and  the 
reference  crop  is  specifically  noted.  Etr  cannot  be  simply 
described  for  all  climate  and  crop  situations.  Relative  leaf 
area  and  the  morphological  and  physiological  aspects  of  the 
crop  canopy  affect  the  energy  exchange  and  the  aerodynamic 
diffusion  process  occurring  above  a  field  (Wright,  1981). 

Two  theoretical  approaches  to  evaporation  from  saturated 
surfaces  are  incorporated  into  the  equations  discussed  here. 
In  the  energy  balance  approach,  evaporation  is  considered  to  be 
a  change  of  state  demanding  a  supply  of  energy  as  heat  of 
vaporization  (Penman,  1956).  The  aerodynamic  approach  assumes 
that  evaporation  is  due  to  turbulent  transfer  of  vapor  by  a 
process  of  eddy  difussion  (Penman,  1948).  Penman  developed  an 
equation  combining  both  of  these  principles  to  predict  "the 
amount  of  water  transpired  in  unit  time  by  a  short  green  crop, 
completely  shading  the  ground,  of  uniform  height  and  never 
short  of  water"  (Penman,  1956). 

Jensen   et   al.   (1971)   described  a  form  of   the   Penman 
equation   which  they  used  in  the  USDA   -  Agricultural   Research 
Service   Irrigation   Scheduling  Program.    The  equation   is   as 
f ol lows : 
E   =  [A/(A+y>]  (Rn  -  G)  +  [y/(A+r)]  (15.36)(1.0  +  0 . 01W) ( es-ed) 


in   which   E   represents  the  daily  potential   evaporative   flux 
from   a   well   watered  reference  crop  like   alfalfa    (Medica go 
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sativa  L.)  with  30  to  46  cm  of  top  growth.  Rn  is  daily  net 
radiation  in  calories  per  square  cm,  G  is  the  daily  soil  heat 
flux  in  cal  per  square  cm,  A  is  the  slope  of  the  saturation 
vapor  pressure  -  temperature  cnrve  (de/dT),  y  is  the 
psychrome trie  constant,  es  is  the  mean  saturation  vapor 
pressure  in  mb  (mean  of  the  saturation  vapor  pressures  at  the 
maximum  and  minimum  daily  air  temperature),  ed  is  the 
saturation  vapor  pressure  at  mean  dew  point  temperature  (or  dew 
point  temperature  near  0800  h)  in  mb,  and  W  is  total  daily  wind 
run  in  miles.  The  parameters  A/fA+y)  and  y/(A+y)  are  mean  air 
temperature   weighting  factors  whose  sum  is  1.0. 


In  their  study  of  peak  water  requirements  for  crops  in 
southern  Idaho,  Wright  and  Jensen  (1972)  realized  that  the 
equation  above  greatly  underestimated  alfalfa  Et .  These  low 
estimates  occurred  primarily  when  a  high  proportion  of  the 
energy  used  for  Et  came  from  the  advection  of  sensible  heat  to 
the  irrigated  area  from  the  surrounding  desert  lands.  They 
therefore  calibrated  the  wind  function  using  2  years  of 
lysimeter  data  from  well  watered,  actively  growing  alfalfa  with 
20  cm  or  more  of  growth.   The  new  wind  function  is  as  follows: 

Wf  =  (0.75  +  0.0185W)   For  W  in  miles/day  at  2  meters. 
This   modification  yields  values  considerably  greater  than   the 
equation  given  previously  when  W  >  50  miles  per  day. 

While  the  Penman  equation  is  favored  by  such  researchers 
as  Wright  and  Jensen  (mentioned  above)  and  Shaw  (1963),  who 
found  it  to  be  the  most  accurate  of  the  methods  he   tested,   it 
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has  some  drawbacks.  Humidity  and  windspeed  data  are  not  widely 
available.  Also,  Hill  et  al.  (1983)  warned  that  judgement 
should  be  used  when  using  any  wind  term  coefficients  if  there 
is  a  different  advective  energy  condition  at  night  than  during 
the  day.  Also,  they  noted  that  field  checks  have  indicated 
that  field  crop  depletion  does  not  follow  the  equation  under 
extremely  windy  conditions.  An  arbitrary  limit  of  100  miles 
per  day  for  W  at  2m  is  imposed  in  some  regions. 


In  addition  to  humidity  and  windspeed,  net  radiation  data 
are  not  widely  available.  Fortunately,  there  are  various 
equations  available  for  estimating  this  input  of  the  Penman 
equation.  Net  radiation  flux  (Rn)  may  be  calculated  from  the 
radiation  balance  equation: 

Rn  =  (1.0  -a)Rs  -  Rln 
in  which  a  is  albedo,  Rs  is  the  global  solar  radiation  flux, 
and  Rln  is  the  longwave  radiation  flux  (de  Jong  et  al.,  1980). 
In  addition  to  equations  for  estimating  solar  radiation  (Rs), 
which  is  readily  available  in  the  areas  in  which  this  research 
was  done,  de  Jong  et  al.  (1980)  presented  two  equations  for 
estimating  Rln.  The  first  was  an  equation  developed  by  Linacre 
(1968),  as  presented  by  de  Jong  et  al.  (1980): 

Rln  =  32  *  10-5  (i  +  4n/N)(100  -  T) 
in   which  T  is  air  temperature  in  degrees  Celsius  and  n  is   the 
number  of  recorded  hours  of  bright  sunshine  in  a  daylength  of  N 
hours.    The   second   equation,   which   they  found   led   to   an 
overestimation   of   Rln  due  to  its  lack  of  a   cloud   adjustment 
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factor,   was  developed  by  Idso  and  Jackson  (1969),  (as  reported 
by  de  Jong  et  al.,  1980): 

Rln  =  a    TK4  (0.261  exp(-7.77  •  10-4  (273  -  Tg)2)). 
In  the  above  equation,  a    is  the  Stefan  -  Boltzmann  constant  and 
T_  is  the  air  temperature  in  degrees  kelvin. 

Wright  (1982)  presented  an  equation  for  net  outgoing 
longwave  radiation  (Kb)  as  an  input  of  the  radiation  balance 
equation  to  calculate  Rn.  The  equation  is: 
Rb  =  (  a(Rs/Rso)  +  b  )Rbo 
in  which  Rs  is  the  incident  solar  radiation,  Rso  is  the  clear 
day  solar  radiation,  Rbo  is  the  net  clear  day  outgoing  longwave 
radiation,  and  the  coefficients  are  dependent  upon  the  ratio  of 
Rs  to  Rso.  An  equation  for  Rbo,  with  the  inputs  of  saturation 
vapor  pressure  at  mean  dewpoint  temperature  (ed)  and  average 
daily  temperature,  is  also  contained   in  Wright's  article. 

Locally  calibrated  Rn  equations  resulted  from  work  done  at 
the  Evapotranspiration  Research  Field  14  km  southwest  of 
Manhattan,  KS.  Rosenthal  et  al.  (1977)  reported  a  relationship 
between  Rn  and  Rs  specifically  for  corn: 


Rn  =  0.861Rs  -  103.92   For  LAI  <  3.0 

Rn  =  0.848Rs  -  144.49   For  LAI  >  3.0 

Rn  -  0.766Rs  -  99.89    For  LAI  <  3.0 

and  after  blister  stage 

where   Rn   and  Rs  are  in  langleys/day  and  LAI  is  the  leaf   area 

index.    Likewise,  such  a  relationship  for  soybean  was  reported 

by  Kanemasu  et  al.  (1976): 


IS 


Rn  =  0.725JU  -  0.86   For  LAI  <  3 
Rn  =  0.805Rs  -  2.332  For  LAI  >  3. 

In  the  above  equations,  Rn  and  Rs  are  in  mm/day. 

The   above   net   radiation  equations   require   either   the 

measurement  or  estimation  of  leaf  area  index.    Various  methods 

of   estimating   soybean  LAI  have  been   proposed.    Wiersma   and 

Bailey   (1975)   suggested  that  total  leaf  area  (TLA)   could   be 

predicted   by  a  linear  equation  having  the  sum  of  the   products 

of   the  lengths  and  widths  of  all  the  terminal  leaflets  as   the 

independent  variable: 

TLA  =  6.532  +  2.045  (  £  LiWi  terminal  leaflets). 

Sivakumar   (1978)  regressed  leaf  dry  weight  and  leaflet   number 

against   measured   leaf   area   per  plant   (LA)   to   derive   the 

prediction  equation: 

LA  =  -286.7  +  80.3(LDW)  +  31.5(LN) 

As  Ogbuehi  and  Brandle  (1981)  pointed  out,   no  independent  test 

of   this   equation  was  reported.    They  too   found   significant 

linear   regressions   of  leaf  area  on  leaf  dry  weight   and   leaf 

number.    While  these  equations  were  accurate  in  predicting   LA 

of   independent   samples   taken   from  the   same   treatment   and 

season,   they  were  inaccurate   in  predicting  LA  of  plants  grown 

in   a  different  treatment  or  season.    Koller   (1972)   reported 

that   the   leaf  area  to  leaf  dry  weight  ratio  changes  over   the 

growing   season  of  the  soybean  plant  and  suggested  that   models 

for  the  prediction  of  soybean  leaf  area  from  leaf  weight  should 

include  time  as  a  variable.    A  different  approach  by   Sinclair 
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(1984)  produced  a  model  of  total  plant  leaf  area  (L)  based  on 
plastochron  index  (PI),  defined  as  an  integer  connt  of  the 
number  of  emerged  leaves  pins  a  decimal  fraction  representing 
the  progress  of  the  emerging  leaf  toward  full  emergence.  The 
relationship  is  as  follows: 

L=A*PI+B»  (exp(C  •  PI(3/2))  -  1) 
in  which  the  coefficients  differ  among  cnltivars. 

Leaf  area  index  development  is  very  dependent  on 
temperature  variations  among  different  environments. 
Goldsworthy  (1975)  compared  the  LAI  vs.  days  after  sowing 
curves  of  a  highland  corn  hybrid  grown  at  an  elevation  of  2,250 
m  above  sea  level  and  a  lowland  variety  grown  at  60  m.  The 
growing  season  at  the  high  elevation  site  was  193  days  while  at 
the  low  elevation  site  it  was  112  days.  Maximum  LAI  was 
attained  at  60  days  after  sowing  in  the  lowland  site  and  at  91 
days  after  sowing  in  the  highland  site.  While  LAI  declined 
rapidly  after  silking  at  the  lowland  site,  at  the  highland  site 
its  near  -  maximum  value  was  maintained  for  an  additional  45 
days  after  silking.  Likewise,  Eik  and  Hanway  (1965)  reported 
that  not  only  did  longer  season  corn  hybrids  produce  larger 
total  leaf  area  per  plant,  but  they  also  tended  to  maintain 
green  leaf  areas  for  longer  periods  than  the  shorter  season  hy- 
brids. 


The  frequent  lack  of  available  data  for  the  Penman 
equation  was  dealt  with  by  Merva  and  Fernandez  (1982).  The 
results     of    their    sensitivity    analysis    showed     that 
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extraterrestrial   radiation  has  the  greatest  influence   of   any 

parameter   on  the  Et  calculation.    Conversely,   a  variation  in 

either   the   dewpoint   temperature  or   wind   velocity   makes   a 

relatively      small     contribution     to     the      predicted 

evapotranspi ra t i on.    They  found  that  a  20%  variation  in  either 

of   these  two  parameters  does  little  to  alter  the  predicted  Et , 

with  wind  contributing  the  least  to  the  variation. 

The   Jensen-liaise   equation  for   estimating   potential   Et 

relies   on  the  following  basic  principles  of  the  energy  balance 

concept:   "When   an  evaporating  surface.   such  as   an   actively 

growing  crop,   is  supplied  with  adequate  water,   the  rate  of  Et 

is   controlled  by  the  available  heat  energy"  (Jensen  and  Haise, 

1963).    Using   data  from  crops  with  adequate  soil  moisture   as 

well   as   evaporating  and  transpiring  surfaces  which   were   not 

limiting  the  vaporization  of  water,   Jensen  and  Haise  developed 

the  following  linear  relationship: 

(Et/Rs)   =  (0.014JT  -  0.37 
P 

In  the  above  equation,  (Et/Rs)  is  the  potential  ratio  of  Et 
to  short  wave  solar  and  sky  radiation  flux  and  T  is  the  mean 
air  temperature  in  degrees  F.  The  equation  can  be  rewritten  to 
estimate  potential  Et,  which  the  authors  define  as  the  Et  which 
can  occur  in  irrigated  fields  located  in  arid  and  semi-arid 
areas.  They  further  explain  that  Etp  does  not  imply  a 
homogeneous  or  unlimited  boundary  area  of  well  watered  actively 
growing  vegetation.  Jensen  et  al.  (1970)  presented  a  modified 
version  of  the  Jensen-liaise  equation  containing  an  air 
temperature  coefficient,  CT#  which  is  constant  for  a  given  area 
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and    derived    from   long-term   mean    maximum    and    minimum 

temperatures   for   tie  month  of  highest  mean   air   temperature. 

Due   to  the  large  changes  in  the  air  temperature-net   radiation 

relationships  at  high  elevations,   one  of  the  constants  used  to 

calculate   CT  is  adjusted  for  elevation.    The  modified  Jensen- 

Haise  equation  takes  the  form: 

Etp  =  C_(T  -  T  )Rs 
T        i 

in  which  T  is  the  mean  daily  air  temperature,  T  is  a  constant 
for  a  given  area  and  is  the  linear  equation  intercept  on  the 
temperature  axis,  and  Rs  is  the  daily  solar  radiation  expressed 
as  equivalent  depth  of  evaporation.  Etp  represents  the  maximum 
Et  that  can  occur  under  given  climatic  conditions  from  a  well 
watered,  aerodynamical ly  rough  crop,  such  as  alfalfa  with  about 
30  cm  to  46  cm   of  top  growth. 

CT  and  T  can  be  determined  by  local  calibration.  They 
can  also  be  calculated  from  the  following  equations  listed 
conveniently  by  Hill  et  al.  (1983): 

cT  =  i/cCj  +c2cH) 

Cx  =  68  -  (3.6) (elev./lOOO) 

C2  =  13 

CH  =  50(e2  -  e,) 

Tx  =  27.5  -  0.25(e2  -  e^  -( el ev . / 1000) 


In   the   above   equati 


and   6j   are   saturation   vapor 


pressures  (mb)  of  water  at  the  long  term  mean  maximum  and 
minimum  temperatures,  respectively,  for  the  warmest  month  of 
the  year  for  the  study  area.    Elev.  is  elevation  in  feet  above 
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sea  level.  Hill  et  al.  suggested  that  caution  be  used  when 
applying  C  and  T  with  full  elevation  correction  for  areas  of 
4500  feet  and  above,  depending  on  the  surroundings.  They 
regressed  correction  factors,  defined  as  the  attainable  Et 
divided  by  the  calculated  Et,  against  elevation,  thus  obtaining 
the  following  linear  equation: 

CFJHE  =  1.653  -  0.1640  E1./1000. 
The  above  equation  relates  a  correction  factor  for  the  modified 
Jensen-Haise  method  to  elevation  in  feet  above  sea  level. 

Hill  et  al.  (1983)  defined  the  Priestley  -  Taylor  equation 
as  essentially  the  Penman  equation  with  both  of  the  wind  term 
coefficients  equal  to  0.   It  takes  the  form: 

Etmax  =  a  [s/(s  +Y)]Rn 
in  which  Etmax  is  daily  maximum  evapot r ansp i ra t i on  during 
predominantly  nonadvective  conditions,  s  is  the  slope  of  the 
saturation  vapor  pressure  curve  at  a  weighted  average 
temperature  3[(Tmax  +  Tmin)/4],  v  is  the  psy chrome tri c 
constant,  and  Rn  is  net  radiation  in  ly/day  (Rosenthal  et  al. 
1977) .  The  coefficient  a  is  dependent  on  climate  and  crop 
type. 

Various  values  for  a  exist  in  the  liturature.  Using  lysi- 
meter,  fluxatron,  and  open  water  data,  Priestley  and  Taylor 
(1972)  found  an  overall  mean  a  of  1.26.  (The  lysimeter  data 
alone  produced  an  o  of  1.32.)  Davies  and  Allen  (1973) 
calculated  values  of  a  which  ranged  between  1.16  and  1.36  and 
had  an  overall  mean  of  1.27.  Mukammal  and  Neman  (1977)  found 
that   near   field  capacity,   a  has  a  value  of  1.29.    De   Bruin 
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(1983)  found  an  a  of  1.3  for  well  watered  surfaces.  Work 
done  at  the  Evapo t r ansp i r a t i on  Research  Field  near  Manhattan 
resulted  in  an  a  of  1.35  for  corn  and  1.45  for  soybean 
(Kanemasu  et  al.  1976:   Rosenthal  et  al.,  1977). 


Heat  Units 

Scientists  have  long  known  that  any  attempt  to  measure 
plant  development  nsing  calendar  days  alone  is  subject  to 
environmental  variation.  As  early  as  1735,  Reaumur  discovered 
that  plant  development  is  more  closely  related  to  the 
temperature  accumulated  to  a  given  stage  than  with  time  alone 
(Neild  and  Seeley,  1977).  While  the  majority  of  researchers 
agree  that  using  "accumulated  heat  units"  or  "growing  degree 
days"  to  predict  phenological  events  is  a  great  improvement 
over  calender  days,  there  is  still  much  controversy  regarding 
the  appropriate  method  of  measuring  heat  units. 

The  oldest  method  was  used  successfully  by  Neild  and 
Seeley  (1977).  The  equation,  GDD  =[(Max  +  Min)/2]  -  50  ,  says 
that  the  growing  degree  days  for  a  given  day  are  equal  to  the 
daily  temperature  average  in  °F  minus  50  °F  (Neild  and  Seeley, 
1977).  The  base  temperature  of  50  °F  is  used  because 
Lehenbauer  (1914)  discovered  that  growth  of  corn  is  extremely 
slow  below  50  °F  (Newman,  1971).  A  linear  regression  of 
Hanway's  numerical  stages  of  development  for  corn  (Ritchie  and 
Hanway,  1982)  against  growing  degree  days  accumulated  from 
planting  to  each  stage  resulted  in  coefficients  of 
determination   (r  )  of  0.98  -  0.99  and  standard  errors  of   0.34 
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to  0.40.    This  compared  favorably  with,  standard  errors  of  0.60 
0.83   when  the  regression  was  done  using  days  from   planting 
(Neild  and  Seeley,  1977)  . 

The  method  described  above  has  been  criticized  on  both 
biological  and  statistical  grounds.  Gilmore  and  Rogers  (1958) 
pointed  out  that  this  method  does  not  represent  the  effective 
heat  units  when  temperatures  fall  below  the  minimum  for  growth 
or  rise  above  the  optimum.  Temperatures  below  50  F  may  not 
further  reduce  growth  and  may  cause  unreasonably  low  GDD 
accumulations.  Relying  on  the  parabolic  nature  of  Lehenbauer's 
growth  curve  for  corn  seedlings,  they  asserted  that 
temperatures  above  the  optimum  retard  growth,  and  therefore 
cause  unreasonably  high  GDD  accumulations.  In  addition,  of  the 
15  methods  of  calculating  heat  units  tested,  they  found  the 
[(Max  +  Min)/2]  -  50  method  to  have  the  highest  coefficient  of 
variation  (CV)  (Gilmore  and  Rogers,  1958).  Likewise,  Mederski 
et  al.  (1973)  found  this  method  to  have  the  highest  CV  of  the  6 
methods  they  tested.   Cross  and  Zuber  (1972)  found  that  when   a 

regression  of  number  of  days  to  pollen  shed  against  accumulated 

2 
thermal  units  was  performed,   this  method  yielded  an  r   of  only 

0.208.    In  contrast,   of  five  models  tested  by  Daughtry  et  al. 

(1984)   ,  including  the  "heat  stress"  method  favored  by  both  the 

teams   of  Gilmore  and  Rogers  and  that  of  Cross  and   Zuber,   the 

daily   average   minus  50  °F  method  produced  the  lowest  CV   when 

used  for  the  period  from  planting  to  silking. 

The  limitations  of  the  standard  GDD  formula  are  dealt  with 
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in  the  modified  growing  degree  day  method  (MGDD)  used  by  the 
National  Weather  Service.  As  far  as  the  calculations  are 
concerned,  any  temperature  below  50  °F  is  set  to  50  °F  and  any 
temperature  above  86  °F  is  set  to  86  °F.  Otherwise,  the 
calculation  is  the  same  as  the  standard  GDD  formula  (Newman, 
1971).  A  similar  method  called  the  "Daily  Adjusted  Average 
System''  had  an  r  of  0.966  when  the  regression  method  of  Cross 
and  Zuber  (as  described  previously)  was  performed.  This  method 
differs  from  the  MGDD  method  in  that  the  base  temperature  is 
not  subtracted  each  day  (Cross  and  Zuber,  1972). 

Another  method,  commonly  called  the  Heat  Stess  Method, 
assumes  that  temperatures  above  the  optimum  retard  growth  and 
should  consequently  result  in  reduction  in  heat  units.  Gilmore 
and  Rogers  (1958)  make  this  correction  by  subtracting  from  the 
daily  mean  the  number  of  degrees  by  which  the  maximum  daily 
temperature  exceeds  the  optimum.  As  in  the  MGDD  method, 
temperatures  below  50  °F  are  considered  to  be  50  °F  and  a  base 
temperature  of  50  F  is  subtracted  from  each  daily  average.  Of 
the  15  methods  tested,  they  found  this  one  to  have  the  lowest 
CV  (Gilmore  and  Rogers,  1958).  Of  the  daily  methods  tested. 
Cross  and  Zuber  (1972)  found  their  daily  heat  stress  method  to 
be  the  best,  but  it  was  only  slightly  better  than  the  daily 
adjusted  average  system.  Their  version  differs  from  the  method 
of  Gilmore  and  Rogers  (1958)  in  that  no  base  temperature  is 
subtracted  and  the  high  temperature  corrections  are  made  by 
subtracting  from  86  °F  the  amount  by  which  the  maximum 
temperature  exceeds  86  °F  and  then  computing  the  average  (Cross 
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and  Zuber,  1972) . 

The  heat  unit  method  commonly  used  in  Canada  is  the 
Ontario  System.  Mederski  et  al.  (1973)  found  this  system  to  be 
the  best  among  those  they  tested  and  describe  it  as  follows: 
Daily  HU  total  =  (day  +  night)/2  in  which  day  =  1.85(max  -  10 
°C)  -  0.026(max  -  10  °C)  and  night  =  min  -  4.4  °C.  A  later 
version  described  by  Coelho  and  Dale  (1980)  was  found  to  be 
less  effective  than  the  HGDD  method. 

Several  researchers  have  attempted  to  devise  a  system 
based  on  Lehenbaner's  growth  curves.  Gilmore  and  Rogers 
(1958)  calculated  "optimum  days"  from  a  freehand  curve  based  on 
Lehenbaner's  data  for  corn  seedlings  grown  at  constant 
temperatures  for  6-hour  periods.  They  fitted  a  scale  to  the 
curve  so  that  the  optimum  temperature  had  a  rating  of  1.00  and 
other  temperatures  had  a  rating  of  a  fraction  until  the  zero 
point  was  reached  when  no  growth  occurred.  They  achieved  a  low 
CV  of  1.24%  by  using  the  mean  of  the  maximum  and  minimum 
temperature  ratings.  Coelho  and  Dale  (1980),  using  a  similar 
technique,  derived  function  of  temperature  equations: 


FT  =  0.027T  -  0.162:  6C  <=  T  <  21C 

FT  =  0.086T  -  1.410:  21C  <=  T  <  28C 

FT  =  1.0:  28C  <=  T  <  32C 

FT  =-0.083T  +  3.67:  32C  <=  T  <  44C 

FT  =  0  for  6C  >  T  =>  44C 

Daughtry  et  al.   (1984)  found  that  differences  in  CV  among 

the   thermal  models   tested  were  very  small   for   planting   to 

physiological  maturity.     However,    for   the    silking    to 

physiological  maturity  interval,   CV  for  calendar  days  was  much 
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smaller  than  CV  for  the  thermal  models.  They  suggested  the  use 
of  a  "mixed"  model  in  which  a  thermal  model  is  used  to  predict 
silking  date  and  then  the  mean  interval  in  days  from  silking  to 
physiological  maturity  is  added  to  predict  date  of 
physiological  maturity.  They  found  that  the  accuracies  of 
these  mixed  models  are  better  than  the  accuracies  of  the 
conventional  thermal  models  for  predicting  physiological 
maturity  of  corn. 

Stapper  and  Arkin  (1980)  presented  a  method  of  calculating 
heat  units  for  corn  which  was  used  in  their  CORNF  growth  and 
development  model.  The  CORNF  growing  degree  day  calculation  is 
identical  to  the  modified  growing  degree  day  method  when  the 
minimum  daily  temperature  is  10  °C  or  above.  However,  when  the 
minimum  daily  temperature  drops  below  10  °C,  a  sine  curve  is 
used  to  approximate  the  diurnal  change  in  temperature  between 
maximum  and  minimum.  The  CORNF  growing  degree  days  are 
multiplied  by  a  daylength  correction  factor,  since  the  actual 
average  temperature  is  higher  than  the  average  of  the  maximum 
and  minimum  when  days  are  longer  than  nights. 

In  comparison  with  corn,  there  is  little  research  in  the 
use  of  thermal  units  for  predicting  soybean  development.  This 
is  most  likely  due  to  the  difficulty  of  accounting  for  the 
photoperiodic  response  of  soybean  in  thermal  models.  Brown 
(1960)  grew  soybean  in  growth  chambers  at  different 
combinations  of  temperature  and  daylength,  finding  that  nearly 
the  same  number  of  night  hours  were  required  to  reach  flowering 


28 


for  all  phot operiods  at  any  one  temperature.  Defining  rate  of 
development  as  the  reciprocal  of  night  hours  from  planting  to 
flowering,  he  plotted  the  rate  of  development  for  each  average 
temperature.  The  curves  thus  derived  were  quadratic 
expressions,  intersecting  the  x  axis  at  50  °F  and  having  an 
optimum  temperature  between  85  °F  and  87  °F  (similar  to  the 
temperature  responses  of  corn). 

Major  et  al.  (1975a)  tested  two  thermal  unit  methods  based 
on  Brown's  quadratic  equation  along  with  9  other  thermal  unit 
methods  similar  to  those  previously  discussed  for  corn. 
Brown's  method,  and  thermal  unit  methods  having  a  base 
temperature  of  10  °C  subtracted  daily,  performed  better  than 
calendar  days  and  methods  with  no  base  temperature.  However, 
no  method  tested  adequately  predicted  post-flowering  events, 
showing   the  need  to  incorporate  daylength  into   a   development 

model . 

In   a  companion  paper,   Major  et  al.   (1975b)  described   a 

mathematical   equation   for   soybean  development   based   on   an 
iterative   regression  analysis  (IRA)  technique.    The  functions 

were   chosen   so  that  the  rate  of  development   was   curvilinear 

with   respect  to  daylength  and  temperature.    The  equation   for 

maturity  (M)  is: 


I         [ax(L  -  ao)  +  a2(L  -  aQ)2]  •  [b^T  -  bQ)  +  bj( 


T  -  b  l"1] 


in  which  L  is  daylength  in  hours  from  sunrise  until  sunset   and 
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T  is  the  average  of  daily  maximum  and  minimum  temperatures  in 
C.  The  a  and  b  coefficients  vary  with  variety  and 
phenological  stage.  S1  and  S-  represent  two  phenological 
stages.  They  found  that  while  the  GDD  method  was  similar  to 
the  IRA  method  for  predicting  the  planting  to  emergence  and 
emergence  to  flowering  periods,  the  IRA  method  was  better  for 
predicting  the  post  flowering  periods  based  on  the  much  smaller 
standard  deviation  of  its  predictions.  This  method  was  used  by 
Curry  et  al.  (1975)  in  SOYMOD  I,  their  model  of  soybean  growth 
and  development. 

Sierra  (1977)  described  a  development  model  for  medium  and 
late  soybean  cultivars  which  was  useful  in  parts  of  Argentina. 
In  order  to  develop  the  model  for  the  emergence  to  flowering 
phase,  a  correlation  and  regression  analysis  was  carried  out 
considering  radiation  summations  as  a  dependent  variable  and 
average  phase  photoperiod  as  the  independent  variable.  This 
regression  revealed  that  the  contribution  towards  flowering  of 
1  langley  received  in  phot oi nduc t i ve  conditions  is  equivalent 
to  progress  towards  flowering  of  2.3S  langleys  received  in 
non-phot oinduct ive  conditions.  Daily  radiation  values  were 
therefore  multiplied  by  2.00  and  used  to  compute  Phot oenerge ti c 
Unit  Summations  (PEUS).  Energetic  -  Phot ot herma 1  Bnit 
Summations  (EPTUS)  were  then  computed  using  the  equation: 

EPTUS  =  PEUS  *  Q-1 
in   which  Q  is  calculated  by  a  complex  equation  having   average 
phase   temperature,   variability  of  average  phase   temperature. 
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and  average  phase  temperature  range  as  variables.  For  the 
flowering  to  ripeness  phase.  Energetic  -  Thermal  Unit 
Summations  were  calculated  as: 

ETUS  =  SRS  *  Z_1 
in   which   SRS   is   the  solar   radiation   summation   and   Z   is 
calculated    from   the   same   variables  as  Q   but   a   different 
equation  is  used. 
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Materials  and  Methods 


Field  Methods. 

Crop  coefficient  curves  based  on  fraction  of  growing 
season  were  developed  for  corn  and  soybean.  These  models  were 
developed  using  data  from  Manhattan  and  Tribune  for  the  years 
1974  through  1982.  The  Manhattan  data  were  collected  at  the 
Evapotranspiration  Research  Field  located  14  km  southwest  of 
Manhattan.  The  Tribune  data  were  collected  at  the  Ross 
Irrigation  Field  near  Tribune.  Specific  information  about  and 
references  for  the  1974  through  1982  corn  and  soybean  plots 
are  contained  in  Tables  1A  and  2A  in  the  appendix.  All 
soybean  cultivars  used  in  this  study  are  indeterminate,  except 
for  the  cultivars  Bay  and  Clark  Determinate.  Two  soils  at 
Manhattan  are  involved  in  this  study.  The  Muir  silt  loam  is  a 
fine-silty,  mixed,  mesic  Cumulic  Haplustoll  (U.S.  Dept.  of 
Agriculture  -  S.C.S.,  1975).  Formed  in  deep  alluvium,  the 
Muir  SiL  is  deep  and  nearly  level  and  has  moderate  permeability 
and  high  available  water  holding  capacity.  The  Eudora  silt  loam 
is  a  coarse  -  silty,  mixed  mesic  Fluventic  Hapludoll  (U.S.D.A.- 
S.C.S,  1975).  Formed  in  coarse  silty  alluvium,  the  Eudora  SiL 
is  deep  and  nearly  level,  with  moderate  permeability  and  high 
available  water  holding  capacity.  The  plots  at  Tribune  were 
located  on  a  Ulysses  SiL,  a  fine  -  silty,  mixed  mesic  Aridic 
Haplnstoll  (Gwin  et  al.,  1974).  The  soil  was  developed  in 
upland  loess.    It  has  a  water  intake  rate  of  about  1.3  cm   per 
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hour   and  0-1%  slopes.    Its  water  holding  capacity  is   high, 
and  its  subsoil  is  calcereous. 

The  crop  coefficient  models  were  tested  using  data 
collected  at  the  Manhattan  site  in  1983  and  1984.  Three  corn 
and  three  soybean  cultivars  were  chosen  for  their  differing 
growing  season  lengths,  since  one  objective  of  the  study  was  to 
evaluate  the  validity  of  basing  crop  coefficient  curves  on 
fraction  of  growing  season.  The  corn  cultivars  used  and  their 
expected  heat  units  needed  for  maturity  (according  to  the 
Cargill  seed  catalogue  for  1983)  were  80S  hybrid,  2040  heat 
units;  872  hybrid,  2410  heat  units;  and  980  hybrid,  2820  heat 
units.  In  1984,  Cargill  hybrid  874  with  2410  heat  units 
replaced  872  in  the  field  plots  (1983-84  Growers  Guide, 
Cargill).  The  soybean  cultivars  used  and  their  respective 
maturity  groups  were  Amsoy  71,  MG  II;  Union,  MG  IV;  and  Bay,  MG 
V. 


In  both  1983  and  1984,  the  study  was  organized  in  a 
randomized  complete  block  design.  The  corn  cultivars  were 
grown  on  the  Muir  SiL  in  1983,  and  the  plot  arrangement  and 
dimensions  are  shown  in  Fig  1.  On  20  April  1983,  224  kg  N/ha 
as  liquid  nitrogen  was  applied  to  the  corn  plots  and  the  corn 
was  planted  on  6  May.  Lasso-Bladez  herbicide  was  applied  the 
evening  after  planting.  On  22  June  1983,  berms  were  built 
around  the  corn  and  soybean  plots  to  prevent  runoff.  The  1983 
corn  plots  were  furrow  irrigated  on  7  July  and  again  on  18 
July.     In  1983,  the  soybean  cultivars  were  grown  on  the  Eudora 
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Fig.  1.   Plot  arrangement  for  corn  in  1983,  showing  plot 
number  (top),  cultivar,  and  block  number  (bottom). 
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SiL  and  the  plot  arrangement  and  dimensions  are  seen  in  Fig.  2. 
Treflan  was  applied  to  this  field  on  4  May  and  the  soybean 
cnltivars  were  planted  on  16  May  1983.  The  soybean  plots  were 
furrow  irrigated  on  21  July   and  again  on  29  July   1983. 

Three  3.66  meter  aluminum  neutron  access  tubes  were 
installed  in  each  corn  plot  in  the  two  center  rows.  Four  such 
tubes  were  installed  in  each  soybean  plot,  also  in  the  two 
center  rows.  Approximately  15  cm  of  each  tube  remained  above 
the  soil  surface,  and  the  tubes  were  kept  plugged  with  rubber 
stoppers  to  keep  out  rain  and  other  material.  Care  was  taken 
not  to  disturb  the  plants  during  installation,  and  any  gaps  in 
the  rows  near  tubes  were  replanted  by  hand  to  maintain  a  good 
stand  in  the  vicinity  of  soil  moisture  measurement. 

Soil  water  content  was  read  from  these  tubes  through  the 
use  of  a  Troxl er  neutron  probe  and  a  Troxler  scaler,  model 
2601.  A  simplified  discussion  of  the  working  of  this  device  is 
as  follows:  The  probe  contains  a  mixture  of  Americium  -  241 
and  Beryllium  as  a  source  of  fast  neutrons.  These  neutrons  are 
emitted  radially  into  the  soil  and  lose  kinetic  energy  as  they 
collide  with  various  atomic  nnclei.  The  average  loss  of  energy 
is  greatest  when  these  fast  neutrons  collide  with  the  hydrogen 
nuclei  of  water,  since  the  two  are  of  nearly  equal  mass.  The 
slowed  neutrons  scatter  randomly  in  the  soil,  and  those 
returning  to  the  probe  are  counted  by  a  detector  cell  filled 
with  BF,  gas.  When  a  slowed  neutron  hits  a  B  nucleus  and  is 
absorbed,  an  alpha  particle  (  a  helium  nucleus)  is  emitted, 
creating  an  electrical  pulse  on  a  charged  wire.    The  number  of 
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Fig.  2.   Plot  arrangement  for  soybean  in  1983,  show- 
ing plot  number  (top),  cultivar,  and  block  number 
( bo  c torn) . 
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pulses  over  a  time  interval  is  counted  by  the  scaler  (Hillel, 
1971).  Using  this  device,  we  measured  soil  water  content  down 
to  3.12  meters  at  approximately  7  to  10  day  intervals.  We 
measured  counts  per  IS  seconds  at  20  consecutive  15.24  cm 
depths  in  the  soil  profile  at  each  tube.  A  15  second  standard 
count  was  taken  within  the  probe  case  of  lead  and  polyethelene 
at  each  tube.  These  rough  data  were  converted  to  mm  of  water 
using  the  calibration  equations  discussed  in  the  calculations 
section  of  this  chapter.  Three  gravimetric  samples  were  taken 
from  each  plot  to  determine  the  water  content  of  the  top  7.62 
cm  of  soil.  These  samples  were  weighed,  dried  in  an  oven  at 
105  C  for  48  hours,  then  reweighed.  Water  content  by  dry 
weight  was  calculated  using  the  formula:  (wet  wt.  -  dry  wt.)/ 
dry  wt.  This  gravimetric  water  content  was  then  multiplied  by 
the  soil  bulk  density,  1.4  g/cm  ,  to  obtain  water  content  by 
v  o  1  urn  e  . 

Plant  samples  were  taken  at  approximately  weekly  intervals 
throughout  the  1983  season  from  at  least  3  rows  into  the  plot 
but  never  from  the  center  six  rows.  One  meter  of  row  was  taken 
from  soybean  and  4  plants  in  a  row  were  taken  from  corn.  Leaf 
area  was  read  using  a  Hayashi  Denko  automatic  area  meter  model 
AAM  -  5.  This  photoelectric  apparatus  measures  the  total  area 
of  objects  by  determining  how  much  the  objects  shade  the 
scanning   light   beam.    Leaf  area  index  (leaf  area  per   ground 

area)  was  calculated  for  soybean  by  dividing  leaf  area  by   7620 

2 
cm  .    This   value  was  obtained  by  multiplying  a  row  spacing  of 
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76.2  cm  by  sample  length  of  100  cm.  For  corn,  leaf  area  was 
divided  by  5604  cm  ,  which  is  a  row  spacing  of  76.2  cm 
mnltiplied  by  a  row  length  of  73.54  cm/4  plants  (calculated 
from  plant  population  taken  at  harvest).  Dry  weights  of 
leaves,  stems,  and  reproductive  parts  were  measured  after  plant 
parts  had  dried  in  an  oven  at  70  °C  for  approximately  1  week. 


Growth  stages  of  the  cultivars  were  monitored  at 
approximately  3  day  intervals.  The  system  outlined  by  Hanway 
and  Thompson  (1971)  was  used  to  determine  growth  stages  in 
soybean.  The  growth  stage  system  decribed  by  Ritchie  and 
Hanway  (1982)  was  used  for  corn. 

Twenty  meters  of  row  were  harvested  from  each  of  the  corn 
plots.  These  samples  were  taken  from  rows  7,  8,  11,  and  12, 
thereby  avoiding  the  rows  containing  neutron  tubes.  Thirty 
meters  of  row  were  harvested  from  each  of  the  soybean  plots 
from  the  two  pairs  of  rows  flanking  the  neutron  tube  rows. 
Plants  were  counted  when  the  samples  were  taken  to  determine 
population.  The  harvest  samples  were  threshed  and  the  kernels 
and  beans  were  weighed.  Then  two  gravimetric  samples  were 
taken  from  each  harvest  sample  and  these  were  treated  in  the 
same  way  as  the  soil  discussed  above.  These  water  contents 
were  used  to  adjust  the  weight  of  soybean  to  13%  moisture  (on  a 
wet  weight  basis)  and  corn  to  15.5%. 

On  4  May  1984  ,  ammonium  phosphate  having  a  grade  of  18- 
46-0  was  applied  to  all  study  areas  at  a  rate  of  308  kg/ha. 
Approximately   224   kg  N/ha  were  knifed  into  the  corn   area   as 
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liquid  nitrogen  solution  (28%)  in  late  April  19  84.  The  1984 
corn  cnltivars  were  grown  on  the  Eudora  SiL  and  the  plot 
arrangement  and  dimensions  are  shown  in  Fig.  3.  The  corn  was 
planted  on  11  May  and  Lasso-Bladex  was  applied  to  this  field 
on  12  Hay.  Three  neutron  access  tubes  were  installed  as 
described  for  the  1983  corn  plots.  Treflan  was  applied  to  the 
soybean  field  on  11  Hay  1984,  and  the  soybean  cultivars  were 
planted  on  18  Hay.  The  1984  soybean  plot  arrangement  and 
dimensions  on  the  Huir  SiL  are  seen  in  Fig.  4.  Three  neutron 
access  tubes  per  plot  were  installed.  On  28  June  1984  the 
berms  were  built  around  the  plots.  The  corn  and  soybean  plots 
were  not  irrigated  in  1984. 

Field  methods  in  1984  were  the  same  as  those  in  1983, 
except  plant  sampling  was  not  done.  Soil  moisture  in  the  1984 
soybean  plots  was  read  with  a  Troxler  neutron  probe  and  scaler 
model  3221.  In  this  case,  a  5  minute  standard  connt  was  read 
every  hour.  Twenty  meters  of  row  were  harvested  from  soybean, 
since  we  had  3  replications  of  each  cnltivar  instead  of  2  as  in 
1983. 
Weather  Data 

Rainfall  amounts  were  measured  at  the  study  sites  with  a 
rain  gauge.  Other  cl ima t ol og i ca 1  data  for  the  10  years 
involved  in  this  study  were  obtained  from  records  at  the 
Weather  Data  Library  of  Kansas  State  University.  Maximum  and 
minimum  temperatures  for  Hanhattan  were  read  with  a  ''max  -  min" 
glass  thermometer.  These  values  for  Tribune  were  taken  from 
" Cl imatol ogical  Data:  Kansas"  published  by  the  NOAA  of  the  D.S. 
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Dept.  of  Commerce. 

Solar  radiation  was  measured  with  an  Epply  PSP  pyranometer 
in  Manhattan  and  was  recorded  at  Tribune  by  a  Licor  pyranometer 
in  a  CE  21  automated  weather  station  (Campbell  Scientific  Inc. 
mi cr ol ogger ) .  Manhattan  windrun  was  measured  with  an 
anemometer  located  at  the  evaporation  station  at  the  Tuttle 
Creek  Resevoir.  This  station  was  maintained  by  the  Army  Corps 
of  Engineers.  Anemometer  data  collected  at  the  Tribune 
experiment  station  was  also  used  in  this  study. 

Previous  to  1981,  0800  h  CST  relative  humidity  (RH)  and 
temperature  for  Manhattan  were  obtained  from  hygrothermograph 
records.  This  information  for  1981  to  the  present  was  obtained 
from  records  of  the  automated  weather  station  ( CR  21)  located 
on  the  research  farm  near  the  K.S.U.  campus.  Tribune  0800  h 
MST  relative  humidity  and  temperature  were  recorded  by  the  CR 
21  located  there . 


£iLl.cul.£t^ons 

The  neutron  probe  readings  were  converted  to  mm  of  water 
by  the  following  procedure.  First,  each  reading  was  divided  by 
the  standard  to  obtain  a  count  ratio  ( CR ) .  The  equation  used 
to  convert  these  count  ratios  obtained  with  probe  and  scaler 
model  2601  to  mm  of  water  is  w  =  (0.45715»CR  -  0.03482)  * 
152.4mm.  The  value  152.4  mm  is  the  thickness  of  each  layer 
whose  water  content  was  measured.  Probe  and  scaler  model  3221 
was  calibrated  by  reading  various  tubes  on  the  same  day  with 
both   3221   and   the   2601   probe   and   scaler.    The   equation 


42 


resulting  from  a  regression  analysis  performed  with  these  data 
is  w  =  (-0.01073  +  0.8033»CR  -  0.25016»CR2)  •  152.4  mm.  The 
3221  probe  and  scaler  was  recalibrated  for  the  1984  soybean 
data  in  the  way  described  above.  The  resulting  equation  is  w  = 
(-0.039055  +  0.90605*CR  -  0.35394*CR2)  •  152.4.  The 
volumetric  water  content  was  multiplied  by  76.2  mm  to  obtain 
the  mm  of  water  in  the  top  layer  of  soil.  The  water  contents 
of  these  layers  were  summed  to  obtain  mm  of  water  down  to  the 
38,  68,  99,  129,  and  160  cm  depths.  For  the  years  1981  through 
1984,  the  water  contents  of  the  190  and  312  cm  depths  were  also 
determined,  as  readings  were  taken  at  deeper  depths  during 
these  years.  The  water  content  sums  were  averaged  for  each 
variety  or  treatment. 

Percent  available  soil  water  was  calculated  for  the  38, 
68,  99,  129,  and  160  cm  layers  by  subtracting  from  measured 
soil  water  content  the  soil  water  content  of  a  given  layer  at  a 
pressure  potential  of  -15  bars,  and  then  dividing  by  available 
water  content  at  field  capacity.  The  -15  bar  water  contents 
for  the  Muir  SiL  were  obtained  from  work  done  by  Brady  (1972), 
and  for  the  Eudora  SiL  these  values  were  obtained  from  work 
done  by  Anderson  et  al.  (1982).  Available  water  content  at 
field  capacity  for  the  soils  was  considered  to  be  the  upper 
limit  of  soil  water  measured  in  the  field  over  the  10-year 
time  period,  minus  the  -15  bar  soil  water  content.  The  values 
thus  obtained  for  the  Muir  SiL  were  80  mm  in  the  38  cm  layer, 
136  mm  in  the  68  cm  layer,  220  mm  in  the  99  cm  layer,  308  mm  in 
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the  129  cm  layer,  and  405  mm  in  the  160  cm  layer.  The  field 
capacity  values  determined  from  the  10  years  of  data  for  the 
Endora  SiL  were  81  mm  in  the  38  cm  layer,  129  mm  in  the  68  cm 
layer,  187  mm  in  the  99  cm  layer,  251  mm  in  the  129  cm  layer, 
and  330  mm  in  the  160  cm  layer.  The  -15  bar  and  field  capacity 
available  water  contents  for  the  Ulysses  SiL  were  derived  from 
work  done  by  Stone  (1982).  These  available  soil  water 
percentages  (IbASW)  were  also  averaged  for  each  variety  or 
t r ea  tment . 

The  Et  rate  for  a  particular  time  period  in  mm/day  was 
calculated  in  the  following  manner.  The  water  content  of  the 
soil  profile  which  was  read  at  the  end  of  the  time  period  was 
subtracted  from  water  content  at  the  start  of  the  time  period, 
and  any  rainfall  amounts  during  the  interval  were  added.  The 
resulting  value  was  then  divided  by  the  number  of  days  in  the 
time  period  to  obtain  mm  of  water  per  day.  These  calculated 
Et  rates  (Etc)  were  averaged  for  each  variety  or  treatment. 
Prior  to  1981,  the  water  contents  in  the  160  cm  layer  were  used 
to  calculate  Et.  For  the  years  1981  through  1984,  water 
contents  in  the  312  cm  layer  were  used  to  determine  Etc. 

Reference  Et  was  calculated  using  the  Penman,  Jensen  - 
Haise,  modified  Jensen  -  Haise,  and  Priestley  -  Taylor 
equations.  The  form  of  the  Penman  equation  used  is  found  in  an 
article  by  Wright  and  Jensen  (1978).  It  estimates  "maximum 
daily  or  potential  Et,  Etp,  for  a  well  -  watered  reference  crop 
of  alfalfa  with  20  cm  or  more  of  top  growth."  The  equation 
takes  the  form: 
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E   ={A/(A+t)  )  (Rn-G)  +  {y/U+Y)}  15.36(0.75  +  0.0115TI)(eo  -e  ) 

z   z 

In   the  above  equation,   E   is  the  estimated  daily   evaporative 

flnx,   Rn   is   net   radiation,   and   G  is   soil   heat   flux   in 

2 
cal/cn  /day.   The  parameter  U  is  wind  run  in  km/day  at  a  height 

of  2  meters;  e    is  the  mean  saturation  vapor  pressure  in  mb  at 

maximum   and  minimum  air  temperature,   and  e   is  the  saturation 

z 

vapor  pressure  based  on  the  0800  h  dewpoint  temperature;  A  is 
the  slope  of  the  saturation  vapor  pressure  -  temperature  curve, 
and  y      is  the  psychrome tr ic  constant  in  mb/°C. 

The  inputs  of  the  Penman  equation  were  calculated  from  our 
weather  data  in  the  following  manner.  Net  radiation  was 
calculated  using  the  locally  calibrated  equations  proposed  for 
corn  by  Rosenthal  et  al.  (1977)  and  for  soybean  by  Kanemasu  et 
al.  (1976).  These  equations  are  listed  in  the  literature 
review.  Aocording  to  Jensen  et  al.  (1971)  "where  day-to-day 
temperatures  do  not  change  greatly  and  day-to-day  radiation  is 
similar,  soil  heat  flux  is  relatively  small  during  the  summer 
months  and  can  be  neglected.''  Since  this  statement  closely 
describes  summer  conditions  in  Kansas,  we  chose  to  disregard 
this  input . 

Our  windrun  was  measured  at  a  height  of  0.61m.  Estimated 
daily  windrun  at  2m,  U2 ,  was  calculated  by  the  following 
equation  recommended  by  Burman  et  al.  (1980): 

U,  =  0   (2/z)0-2 

I  Z 

in   which  z  is  the  elevation  of  the  wind  measurement  and  U    is 

z 

measured   windrun   at   a  height  of   z   meters.    The   foil  owing 
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expression  was  used  to  approximate  A  ,   the  slope  of  the   vapor 
pressure  -  temperature  curve  in  mb/°C: 

A  =  2.00(0.00738T  +  0.8072)7  -  0.00116  (Burman  et  al.,  1980) 
in  which  T  is  mean  daily  temperature  in  °C.  The  following  set 
of  equations  was  used  to  calculate  y  ,  the  psychrometer 
constant  in  mb/°C: 

y  =  0.386P/L 

P  =  1013  -  0.1055E 

L  =  595  -  0.51T   (Burman  et   al.,  1980). 
In  the  above  equations,   P  is  station  barometric  pressure  in  mb 
corrected   for  elevation,   L  is  the  latent  heat  of  vaporization 
in   cal/g,   E   is   sea   level  elevation  in   meters   (314   m   at 
Manhattan  and  1067  m  at  Tribune),    and  T  is  temperature  in   °C 

(we   used   the  average  daily  temperature).    The  input  e°    was 

z 

calculated  by  taking  the  average  of  the  saturation  vapor 
pressure  at  the  maximum  daily  temperature  and  at  the  minimum 
daily  temperature.  These  saturation  vapor  pressures  were 
determined  by  the  following  equation  listed  by  Burman  et  al . 
(1980) : 

es  "=•  33.8639H0.00738T  +  0.8072)8  -  0  .  000019  I  1 .  8T  +  48  I  + 
0.001316] 

in  which  es  is  in  mb  and  T  is  in  °C.  Saturation  vapor  pressure 
at  the  0800  h  dew  point  temperature,  e  ,  was  determined  from 
0800  h  relative  humidity  (RH)  and  temperature  data.  The  0800  h 
temperature  was  used  in  the  above  equation  to  calculate 
saturation  vapor  pressure,  and  this  value  was  then  multiplied 
by   RH   (expressed   as  a  fraction)  to   obtain   e  .    The   daily 
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average  temperature  was  used  to  calculate  the  daily  value  of  L, 

latent  heat  of  vaporization  in  cal/g.    This  value  was   divided 

2 
by   10  to  obtain  L  in  cal/(cn  *mm) ,   and  the  daily   evaporative 

fluxes   calculated  with  all  of  the  reference  Et  equations   were 

divided  by  this  value  to  obtain  Etr  in  mm/day. 

Two   versions   of  the  Jensen  -  Haise  equation  were   used   to 

calculate  Etr.   The  standard  Jensen  -  Haise  equation, 

Etp  =  (0.014T  -  .37)Rs   (Jensen  and  Haise  1963) 

relates  reference  Et  to  solar  radiation  in  equivalent  depth   of 

2 
evaporation   (in    our   case,   cal/cm  /day).    T  is  daily   mean 

temperature   in    F.    The   modified   Jensen   -  Haise   equation 

(Jensen  et  al. ,  1970)  was  also  used  to  calculate  Etr: 

Etr  =  CT(T  -  Tx)Rs    (Burman  et   al.,  1980) 

In  the  above  equation,  T  is  daily  average  temperature  in  °C  and 

2 
Rs  is  solar  radiation  in  cal/cm  /day.    The  constants  C   and  T 

were  calculated  by  the  following  series  of  equations  listed   by 

Burman  et   al.  (1980): 

CT  =  1/tCj  +  7.3CH) 

CH  -  50  mb/(e2  +  e±) 

C1    =  38  -  2E/305 

Tx  =  -2.5  -  0.14(e2  -  e^  -  E/550 
In  the  above  equations,  e,  is  the  saturation  vapor  pressure  in 
mb  at  the  mean  monthly  maximum  air  temperature  of  the  warmest 
month  of  the  year  (33.2  °C  for  Manhattan  and  33.7  °C  for 
Tribune  in  July),  and  e.  is  the  saturation  vapor  pressure  in  mb 
at   the   mean   monthly  minimum  air  temperature  of   the   warmest 
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month   of   tie   year  (20.0  °C  for  Manhattan   and   16.5   °C   for 

Tribune).    E   is   the   site   elevation  in   meters   (314   m   at 

Manhattan   and  1067  m  at  Tribune).    The  modified   Jensen-Haise 

equation   constants   for  Manhattan  were  0.020315  for  C_   and 

6.9177  for  T  .   For  Tribune,  these  constants  were  0.0238712  for 

C_  and  -9.13333  for  T  . 
1  i 

The   Priestley   -  Taylor   equation   used   is   found   in   a 

publication  by  Kanemasu  et  al .  (1976): 

ETmai  =  o  [»/(»+•»}]  Rn 

In   the   above   equation,    ETmax   is   the   " energy- 1 im i te d   Et 

occurring   from   a   well-watered   surface   during   nonadvective 

conditions''   (Kanemasu  et  al.,   1976).    Rn  is  net  radiation  in 

2 
cal/cm  /day   and   was  calculated  as  described   for   the   Penman 

equation.    The   constants   s   and   y   are   the   slope   of   the 

saturation     vapor    pressure-temperature    curve    and     the 

psychrometric  constant  in  mb/°C  at  mean  daily  temperature.    We 

used  an  a  of  1.32  for  both  corn  and  soybean. 

The   time   scale  used  in  the  crop   coefficient,   leaf   area 

index,   and   depletion    curves   presented   in   this   paper   is 

fraction   of   growing  season.    Fraction  of  growing   season   is 

defined   as   the  heat  units  accumulated  from  emergence   to   the 

time   in   question,   divided  by  the  total  number  of  heat   units 

accumulated  from  emergence  to  physiological  maturity. 

In   order   to   choose   a  heat   unit   method   to   calculate 

fraction   of  growing  season  for  corn,   5  heat  unit  methods,   as 

well  as  days  after  emergence,   were  used  to  calculate   fraction 

of   growing  season  to  silking  for  each  year  -  cultivar  data  set 
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(see  Table  3A  in  the  appendix).  This  analysis  is  discussed  in 
the  results  and  discussion  section  of  this  paper.  The  methods 
used  were  the  growing  degree  day  method  (ODD),  the  modified 
growing  degree  day  method  (HGDD) ,  the  Gilmore  and  Rogers  (1958) 
heat  stress  method,  the  Cross  and  Zuber  (1972)  heat  stress 
method,  and  the  CORNF  method.  In  the  following  equations,  Tmax 
is  the  daily  maximum  temperature  in  C  and  Tmin  is  the  daily 
minimum  temperature  in  C.  Tmax  is  30  C  if  Tmax  >  30  C,  and 
Tmin  is  10°C  if  Tmin  <  10°C.  Tavg  =  (Tmax  +  Tmin)/2. 
GDD  =  [(Tmax  +  Tmin)/2]  -  10. 

MGDD  =  [(Tmax*  +  Tmin*)/2]  -  10. 
Heat  Stress  (GE)  =  {[(Tmax  +  Tmin*)/2  ]  -  (Tmax  -  30))  -  10   if 
Tmax  >  30  ,  and   [(Tmax  +  Tmin  )/2]  -  10     if  Tmax  <  30. 
Heat   Stress   ( CZ )  =  {[30  -  (Tmax  -  30)]  +  Tmin*)/2   if  Tmax   > 
30,  and   (Tmax  +  Tmin  )/2    if  Tmax  <  30. 

CORNF  -  1/n  {AMP  •  COSINE(ZETA)  +  (Tavg  -  10)(n/2  -  ZETA) )     if 
Tmin  <  10,  and   {(Tmax   +  Tmin)/2)  -  10    if  Tmin  =>  10. 
AMP  =  Tmax  -  Tavg. 

ZETA  =  ARCSINE[(10  -  TavgWAMP]. 
The   CORNF   heat   units  were  multiplied  by  a   daylength   factor 
(HUDAYL)  calculated  by: 

HUDAYL  =  1  -  (14.2  -  DAYLN)  •  0.10 
in  which  DAYLN  is  daylength  in  hours  and  14.2  is  the  average  of 
the   average   corn  growing  season  daylengths  at   Manhattan   and 
Tribune.    The   CORNF   model   is  described  in   more   detail   by 
Stapper   and   Arrin  (1980).    We  chose  the  growing   degree   day 
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method  (GDD)  as  the  basis  for  fraction  of  growing  season  for 
corn.  In  accumuul a t i ng  GDD  to  calculate  fraction  of  growing 
season,  any  negative  daily  value  was  considered  to  be  0. 

The  selected  method  of  calculating  phot othermal  units 
(PTC)  for  soybean  was  described  by  Major  et  al.  (1975b).  The 
basic  equation  is: 


M   =   X    [a1(L-aQ)  +  a2(L-aQ)2]  •  [b1(T-bQ)  +  bjfT-bJ2] 

Sl 
in  which  M  is  accumulated  PTU,  T  is  the  average  of  the  daily 
maximum  and  minimum  temperatures  in  °C  and  L  is  daylength  in 
hours  from  sunrise  to  sunset.  Major  et  al  .  (1975b)  presented 
the  daylength  and  temperature  coefficients  for  10  cultivars, 
two  in  each  of  the  maturity  groups  I  through  V.  One  set  of 
coefficients  was  used  from  emergence  to  flowering  and  another 
set  was  used  from  flowering  to  physiological  maturity.  Since 
the  present  study  involved  cultivars  from  maturity  groups  II 
through  V,  we  used  the  average  of  the  coefficients  of  the  two 
cultivars  in  each  maturity  group.  We  found  that  the 
recommended  coefficients  for  flowering  to  physiological 
maturity  resulted  in  an  unusually  high  fraction  of  growing 
season  for  emergence  to  beginning  bloom  for  an  MG  IV  cultivar 
in  a  hot  season  (see  Table  4A  in  appendix).  We  decided  to  use 
only  MG  II  coefficients  for  flowering  to  physiological  maturity 
for  all  cultivars.  A  comparison  of  the  performance  of  the 
specific   F-PM  coefficients  and  the  MG  II  F-PM  coefficients   is 
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Table      1.  Daylength      and      temperature      coefficients      used      in 

photothermal    unit    equation   for    soybean. 


Daylength  Temperature 


Phase 
and 

™ a1 a2 bo  b1  b2 

E-F 

II  9.02  0.02503         -0.003095            3.00           0.03952         0.0 

III  9.24  0.02695        -0.003641           3.61           0.04196        0.0 

IV  8.89  0.02598         -0.003472            2.72            0.04052         0.0 

V  17.84  -0.01351            0.0                         6.92            0.02431         0.0 

F-PM 

All    18.07  -0.01932     0.0         11.46     0.03556   -0.0 

01391 
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shown  in  the  results  and  discnssion  section  of  this  paper.  The 
coefficients  which  we  used  in  the  above  equation  are  listed  in 
Table  1.   Daylength  in  hours  was  calculated  from  the  equation: 

L  -  12  +  A  •  COSINEf  0.0161  *  (D  -  172)] 
in  which  A  is  the  number  of  hours  that  daylength  on  the  summer 
solstice  exceeds  daylength  at  the  equinox,  and  D  is  the  day 
number  from  1  January.  In  this  equation,  12  represents  the 
daylength  at  the  equinox,  172  adjusts  the  day  number  to  the 
summer  solstice,  and  0.0161  converts  the  day  number  to  radians 
(Major  et  al.,  1975b).  In  accumulating  PTU  to  calculate 
fraction  of  growing  season,  any  negative  daily  value  was  set  to 
0. 
Model  Devel.025.ejt 

Curves  of  maximum  depth  of  depletion  vs.  fraction  of 
growing  season  were  derived  for  corn  and  soybean.  This  was 
done  by  first  subtracting  total  mm  of  water  in  each  15.24  cm 
layer  (below  the  top  7.62  cm)  on  the  last  day  from  mm  of  water 
in  the  same  layer  on  the  first  day  of  a  time  period.  Depletion 
in  mm  was  calculated  in  this  way  for  each  layer  down  to  the 
deepest  reading  depth  (160  cm  before  1981,  312  cm  from  1981 
onwards).  Maximum  depth  of  depletion  was  considered  to  be  the 
depth  after  which  depletion  dropped  to  a  substantially  lower 
amount.  Maximum  depth  of  depletion  was  determined  only  for 
time  periods  in  which  rainfall  averaged  less  than  1  mm/day. 
All  such  data  from  1974  through  1984  were  used  in  these  models. 
Fraction  of  growing  season  using  the  growing  degree  day  method 
for   corn   and   the  photothermal  unit  method   for   soybean   was 
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calculated  for  the  first  and  last  days  of  a  time  period.  These 
two  numbers  were  averaged,  and  maximum  depth  of  depletion  was 
regressed  against  this  mean  fraction  of  growing  season.  The 
STEPWISE  procedure  of  the  SAS  statistical  analysis  computer 
system  was  used  to  choose  the  best  regression  equations  (SAS 
Istitute  Inc.,  1982).  The  maximum  R  improvement  technique 
(MAIR)  was  used,  and  the  independent  variables  were  fraction  of 
growing  season  raised  to  the  first  through  the  fifth  powers.  I 
considered  the  best  model  to  be  that  with  the  highest  R2  as 
well  as  parameters  which  were  significant  at  the  10%  level. 
This  procedure  was  also  used  for  choosing  the  LAI  and  crop 
coefficient   models.    Occasionally,   a  more  complicated   model 

with  significant  parameters  was  disregarded  if  its   improvement 

„2 
in  R   was  very  small. 


Since,  as  discussed  earlier,  I  calculated  available  soil 
water  for  the  38,  68,  99,  129,  and  the  160  cm  layers,  I  solved 
these  depletion  equations  for  the  fractions  of  growing  season 
at  which  the  equations  equaled  these  depths  of  maximum 
depletion.  These  values  served  as  estimates  of  effective 
rooting  zone  throughout  the  growing  season.  The  ASW  values 
which  had  been  calculated  for  these  layers  on  each  reading  day 
were  used  in  the  crop  coefficient  curves  described  later. 

Since  the  net  radiation  equations  in  this  model  require  an 
estimate  of  when  LAI  reaches  3  and  then  drops  below  3,  we  de- 
veloped LAI  models  using  the  1981  and  1982  data.  Each 
measured   LAI   value  was  divided  by  the  maximum  LAI  reached   by 
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that  particular  ooltivar  in  the  growing  season  of  the 
measurement.  These  values  of  "relative  leaf  area  index"  were 
regressed  against  fraction  of  growing  season  and  a  regression 
equation  was  chosen  using  STEPWISE  ( SAS  Institute  Inc.,  1982) 
as  described  for  maximum  depth  of  depletion. 

These  curves  of  relative  LAI  vs.  fraction  of  growing 
season  were  then  used  to  predict  relative  LAI  for  the  1983 
data.  The  predicted  relative  LAI  values  were  regressed  against 
the  measured  values  using  a  linear  model  in  the  REG  procedure 
in  SAS  (SAS  Institute  Inc.,  1982).  An  F  test  was  used  to  test 
if  the  slopes  of  the  predicted  vs.  measured  lines  equaled  1  and 
if  the  intercepts  equaled  0.  This  was  done  using  the  TEST 
procedure  in  REG.  The  corn  and  soybean  LAI  models  were  solved 
for  the  fractions  of  growing  season  at  which  relative  LAI  is 
equal  to  0.6,  producing  an  LAI  of  3  if  maximum  LAI  is  assumed 
to  be  5.  For  corn,  if  fraction  of  growing  season  was  between 
0  and  0.287,  the  first  Rn  equation  for  corn  listed  in  the 
literature  review  was  used  in  the  Penman  and  Priestley  -  Taylor 
equations.  If  fraction  of  growing  season  was  betwen  0.287  and 
0.831,  the  second  equation  was  used  and  the  third  was  used  for 
fractions  of  growing  season  greater  than  0.831.  For  soybean, 
the  first  equation  was  used  when  fraction  of  growing  season  was 
less  than  0.521,  the  second  when  fraction  of  growing  season  was 
between  0.521  and  0.876,  and  the  first  equation  was  used  again 
when  fraction  of  growing  season  was  greater  than  0.876. 

The  crop  coefficients  for  corn  and  soybean  were  developed 
by   calculating  the  ratios  of  actual  Et  (Etc)  to   reference   Et 
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(Etr),  adjusting  these  ratios  for  percent  available  soil  water, 
and  regressing  these  Kcb's  against  fraction  of  growing  season. 
The  data  used  in  these  basal  crop  coefficient  curves  were  all 
neutron  probe  reading  time  periods  in  which  rainfall  was  less 
than  3  mm  for  the  yean  1974  through  1982.  Actual  Et  (mm/day) 
was  calulated  for  a  given  time  period  in  the  way  described  in 
the  calculations  section.  Reference  Et  for  each  day  was 
calculated  using  the  Penman,  Priestley  -  Taylor,  Jensen-liaise, 
and  modified  Jensen-Baise  equations.  Reference  Et  in  mm/day 
(Etr)  for  a  time  period  was  determined  by  summing  the  daily 
values  from  the  first  day  of  a  time  period  through  the  last 
day,  and  then  dividing  this  value  by  the  number  of  days  in  the 
time   period.    Actual  Et  (Etc)  in  mm/day  was  divided  by  Etr  in 

mm/day   to   obtain  the  crop  coefficient  (Kc)  for  a   given   time 

period. 

These   Kc  values  were  adjusted  for  percent  available   soil 

water  since,   as  described  in  the  literature  review,   available 

water   within  the  effective  rooting  zone  limits   transpiration. 

These    basal   crop   coefficients   (Kcb)   were   calculated    by 

rearranging  the  equation: 

Kc  =  Kcb*Ka  +  Ks    (Burman  et  al . ,  1980) 

to  read: 

Kcb  -  Kc/Ka. 

The   rainfall  correction  term  was  disregarded  since  only   times 

of  no  or  very  little  rainfall  were  included  in  the  model.    The 

ASW  correction  Ka  was  calculated  by: 
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Ka  =  [ln(%ASW  +  1 ) ] / [ 1 n ( 101 ) ] 
in   which   %AStf  is  the  percent  of  available  soil  water   in   the 
effective  rooting  zone  and  In  is  the  natnral  log  function.   The 
%ASW   nsed  was  the  average  of  the  %ASW  values  for  the  first  and 
last  days  of  a  time  period. 

Fraction  of  growing  season  using  the  methods  previously 
described  was  calculated  for  the  first  and  last  days  of  a  time 
period  and  these  two  values  were  averaged.  The  STEPWISE 
procedure  was  performed  on  the  crop  coefficient  curve  data  as 
described  for  the  depletion  curves.  The  MAXR  option  was  again 
used  with  Kcb  serving  as  the  dependent  variable  and  the  first 
through  fifth  powers  of  average  fraction  of  growing  season 
serving  as  the  independent  variables  (SAS  Institute  Inc., 
1982).  The  models  chosen  were  those  which  had  the  highest  R^ 
values  as  well  as  parameters  significant  at  the  10%  level.  In 
the  case  of  soybean,  several  such  models  were  plotted  for  each 
Etr  equation  used,  and  the  most  reasonable  ones  were  chosen. 

The  equations  of  these  derived  Kcb  curves  were  used  to 
predict  Et  for  neutron  probe  reading  time  periods  in  1983  and 
1984.  These  predicted  values  were  then  compared  with  the 
measured  values.  The  average  fraction  of  growing  season  for  a 
time  period  was  the  first  input  used  in  the  basal  crop 
coefficient  curve  equation  to  yield  the  Kcb.  Average  %ASW  for 
the  time  period  was  used  to  calculate  a  Ka  value  as  explained 
above.  The  Kcb  was  then  multiplied  by  Ka  to  obtain  Kc .  This  Kc 
was  then  multiplied  by  Etr  in  mm/day.  This  Etr  was  calculated 
by   the  equation  used  to  derive  the  Kcb  curve  being  used  in  the 
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prediction. 

Since  Et  for  all  tine  periods  in  1983  and  1984  were 
predicted  regardless  of  rainfall,  the  Ks  factor  proposed  by 
Jensen  et  al.  (1971)  was  added  to  the  predicted  Et  if  rain  did 
occur.  Ks  was  approximated  by  (0.9  -  Kcb)0.8,  (0.9  -  Kcb)0.5, 
and  (0.9  -  Kcb)0.3,  for  the  first,  second,  and  third  days  after 
a  rainfall,  respectively  (Burman  et  al.,  1980).  We  considered 
Kcb  for  the  rainfall  corrections  to  be  the  average  Kcb  of  the 
time  period  obtained,  as  described  above,  by  using  average 
fraction  of  growing  season  in  the  crop  coefficient  cnrve 
equation.  The  three  Ks  values  were  totaled  each  time  a 
rainfall  event  occurred,  the  total  amount  of  correction  never 
exceeding  the  actual  amount  of  rain.  Any  leftover  correction 
value  from  one  time  period  was  added  onto  the  next  time 
period.  Before  the  total  of  the  Ks  values  for  a  time  period 
was  added  onto  the  Et  predicted  by  Kc*Etr  (mm/day),  it  was 
divided  by  the  number  of  days  in  the  time  period.  As  described 
for  the  prediction  of  LAI,  the  predicted  Et  values  were 
compared  with  the  measured  values  using  the  REG  procedure  in 
SAS.  Linear  regression  was  performed  and  the  TEST  option  was 
used  to  test  equality  to  1  of  the  slopes  of  the  predicted  vs. 
measured  Et  lines.  TEST  was  also  used  to  test  equality  of  the 
intercepts  to  0  (SAS  Institute  Inc.,  1982). 
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Results  and  Discussion 


The  leaf  area  index,  depletion  depth,  and  crop  coefficient 
curves  presented  in  this  section  are  based  on  fraction  of 
growing  season  (FGS)  as  their  time  scale.  For  corn,  this 
fraction  was  computed  using  the  growing  degree  day  method  (GDD) 
described  in  the  last  chapter.  This  method  was  chosen  by 
analyzing  fraction  of  growing  season  to  silking,  calculated  by 
using  six  methods,  for  18  y ear-cul t i var  data  sets  (see  Table 
3A) .  I  used  the  coefficient  of  variation  (CV)  in  evaluating 
various  heat  unit  methods.  I  intended  to  combine  data  sets  of 
different  years,  locations,  and  cultivars  into  FGS-based 
models,  and  then  use  fraction  of  growing  season  as  an  input 
into  the  models  in  order  to  make  predictions  about  independent 
test  data.  Consistency  of  results  was  therefore  a  criterion 
in  choosing  heat  unit  methods.  CV,  estimated  by  sample 
standard  deviation  divided  by  the  sample  mean,  is  a  relatively 
stable  measure  of  variation  because  it  is  independent  of  the 
sample  mean.  It  is  therefore  a  suitable  measure  of  consistency 
of  results  (Snedecor  and  Cochran,  1980). 

Tables  2  and  3  show  this  analysis  for  corn,  performed  with 
and  without  the  1983  and  1984  data  used  to  test  the  crop 
coefficient  curves.  The  CV  of  the  GDD  method  is  only  slightly 
higher  than  those  of  the  other  heat  unit  methods.  Considering 
also  the  simplicity  of  the  GDD  method,  I  found  it  to  be  the 
most  attractive  of  the  methods  tested. 
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Table  2.  Fraction  of  growing  season  to  silking  in  corn, 
calculated  by  using  five  heat  unit  methods  and  days 
after  emergence  (DAE).  The  analysis  did  not  include 
the    1983    and    1981    data.       N=12. 


Method 

Fraction 

Std.    dev. 

cv 

(*) 

GDD 

0.502 

0.039 

7.85 

MGDD 

0.508 

0.034 

6.66 

Heat    stress    (GR) 

0.516 

0.032 

6.26 

Heat    stress    (CZ) 

0.522 

0.032 

6.08 

CORNF 

0.531 

0.038 

7.10 

DAE 

0.528 

0.037 

6.  19 
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Table  3.  Fraction  of  growing  season  to  silking  in  corn, 
calculated  by  five  heat  unit  methods  and  days  after 
emergence  (DAE).  The  analysis  included  the  1983  and 
1984    data.  N= 1 8. 


Method 


Fraction 


Std.     dev. 


CV 

(?) 


GDD 

0.500 

0.0  42 

MGDD 

0.515 

0.037 

Heat    stress    (GR) 

0.532 

0.0  40 

Heat    stress    (GR) 

0.539 

0.041 

CORNF 

0.535 

0.0  40 

DAE 

0.548 

0.0  42 

8.30 
7.21 
7.60 
7.53 
7.47 
7.72 
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A  similar  analysis  was  made  of  fraction  of  growing  season 
to  beginning  bloom  calculated  for  22  soybean  cul t i var-y ear  data 
sets  (see  Table  4A) .  In. the  first  method  used  to  calculate 
FGS,  coefficients  specific  to  each  maturity  group  were  used  in 
the  phot othermal  unit  equation  described  in  the  Materials  and 
Methods  section.  The  second  method  involved  using  specific 
coefficients  for  the  emergence  to  flowering  period  and  MGII 
coefficients  only  for  the  flowering  to  physiological  maturity 
period  (F-PM).  The  third  method  used  was  days  after  emergence. 
Tables  4  and  5  show  this  analysis  performed  with  and  without 
the  1983  and  1984  test  data.  Without  the  1983  and  1984  data, 
the  CV  values  are  very  close.  However,  the  addition  of  the  six 
extra  data  sets  considerably  raises  the  CV  for  both  the 
"specific  coefficient"  method  and  the  days  after  emergence 
method.  Since  the  CV  of  the  method  using  only  MGII  F-PM 
coefficients  is  essentially  unchanged  by  adding  the  test  data 
to  the  analysis,  we  considered  this  to  be  a  more  consistent 
method  and  used  it  to  calculate  FGS  for  our  soybean  curves. 

The  curves  of  relative  leaf  area  index  vs.  fraction  of 
growing  season  for  corn  and  soybean  are  shown  in  Figs.  5  and  6. 
Relative  leaf  area  index  is  defined  as  measured  LAI  divided  by 
the  maximum  LAI  achieved  by  the  particular  cultivar  in  the 
growing  season  of  the  measurement.  These  curves  were  developed 
by  using  the  1981  and  1982  data  from  Manhattan  and  Tribune 
described  by  Hattendorf  (1982)  and  Redelfs  (1983)  (see  Table 
1A) .  These  curves  were  then  used  to  predict  corn  and  soybean 
relative   LAI   for   the   1983   data.     Linear   regressions   of 
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Table  4.  Fraction  of  growing  season  to  beginning  bloom 
in  soybean,  calculated  by  using  the  maturity  group- 
specific  coefficients  and  MG  II  coefficients  for 
flowering  to  physiological  maturity  in  the  PTU 
equation,  as  well  as  using  days  after  emergence  (DAE). 
The  analysis  did  not  include  the  1983  and  1984  data. 
N=  1 6  . 


Method 


Fraction 


Std.     dev. 


CV 
(») 


Specific 

coefficients 

F-PM 

0.532 

0.050 

MGII 

coefficients 

F-PM 

0.45  9 

0.0  42 

DAE 

0.332 

0.028 

9.38 

9.05 
8.48 
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Table  5.  Fraction  of  growing  season  to  beginning  bloom 
in  soybean,  calculated  by  using  the  maturity  group- 
specific  coefficients  and  MG  II  coefficients  for 
flowering  to  physiological  maturity  in  the  PTU 
equation,  as  well  as  using  days  after  emergence  (DAE). 
The    analysis   included   the    1983    and    1981    data.    N=22. 


Method 


Fraction  Std.     dev. 


CV 
(?) 


Specific 

coefficients 

F-PM 

0.533 

0.056 

MGII 

coefficients 

F-PM 

0.467 

0.0  40 

DAE 

0.3  49 

0.048 

10.56 

8.61 
13.74 
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Fig.  5.  Regression  curve  of  corn  relative  leaf  area  index 

vs.  fraction  of  growing  season  based  on  growing  degree 

days.   The  curve  was  developed  by  using  the  1981  and  1982 
data  . 
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Fig.  6.  Regression  curve  of  soybean  relative  leaf  area 
index  vs.  fraction  of  growing  season  based  on  photo- 
thermal  units.  The  curve  was  developed  by  using  the 
1981    and    1982    data. 
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predicted  relative  LAI  for  corn  and  soybean,  calculated  by 
using  the  chosen  prediction  equations,  on  measured  relative  LAI 
are  shown  in  Figs.  7  and  8.  My  criteria  for  successful 
prediction  equations  were  slopes  equal  or  nearly  equal  to  1  and 
intercepts  equal  or  nearly  equal  to  0  for  linear  regression 
lines  of  predicted  vs.  actual  values. 

The  TEST  procedure  within  the  REG  procedure  in  the  SAS  stat- 
istical analysis  system  (SAS  Institute  Inc.,  1982)  was  used  to 
test  these  slopes  and  intercepts.  For  the  corn  GDD  method,  the 
null  hypothesis  that  the  slope  of  the  predicted  vs.  measured 
relative  LAI  line  is  equal  to  1  yielded  an  F  value  of  1.034, 
significant  at  the  0.32  level.  The  null  hypothesis  that  the 
intercept  of  this  line  is  equal  to  0  resulted  in  an  F  value  of 
1.40,  significant  at  the  0.24  level.  Since  these  F  values  are 
smaller  than  F  significant  at  the  0.05  level,  the  hypotheses 
that  the  slope  is  equal  to  1  and  that  the  intercept  is  equal  to 
0  can  not  be  rejected.  I  therefore  concluded  that  this  LAI 
curve  for  corn  is  suitable  for  estimating  the  fractions  of 
growing  season  at  which  corn  LAI  is  equal  to  3.0.  This 
information  was  used  in  net  radiation  equations  as  described  in 
the  Materials  and  Methods  section. 

The  soybean  relative  LAI  vs.  FG S  curve  also  proved  to  be 
successful  in  predicting  the  1983  relative  LAI.  The  TEST 
procedure  calculated  an  F  value  of  1.56,  sigificant  at  the  0.22 
level  for  the  null  hypothesis  that  the  slope  of  the  predicted 
vs.   actual  line  is  equal  to  1.    For  the  null  hypothesis   that 
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Fig.  7.  Regression  line  of  predicted  relative  leaf  area 
index  for  corn  vs.  measured  relative  leaf  area  index 
for  the  1983  data. 
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Fig.  8.  Regression  line  of  predicted  relative  leaf  area 
index  for  soybean  vs.  measured  relative  leaf  area  index 
for  the  1983  data. 
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the  intercept  equals  0,  the  procedure  produced  an  F  value  of 
0.01,  significant  at  the  0.94  level. 

I  also  regressed  the  1981  and  1982  relative  LAI  data 
against  days  after  emergence  (DAE)  and  fraction  of  growing 
season  based  on  DAE  and  modified  growing  degree  days  for  corn. 
The  equations  of  these  curves  for  corn  are  shown  in  Table  6, 
and  in  Table  7  for  soybean.  These  equations  were  then  used  to 
predict  the  1983  relative  LAI,  and  the  TEST  procedure  was  used 
to  evaluate  the  predicted  vs.  actual  relative  LAI  regression 
lines.  Table  8  shows  the  predicted  vs.  actual  corn  relative 
LAI  regression  equations  with  their  F  statistics.  Table  9 
shows  this  information  for  soybean.  The  slope  of  the  predicted 
vs.  actual  corn  relative  LAI  line  using  the  modified  growing 
degree  day  method  is  significantly  different  from  1.  The  F 
values  of  the  lines  resulting  from  prediction  equations  for 
corn  based  on  days  after  emergence  and  fraction  of  growing 
season  using  DAE  have  F  statistics  which  are  highly 
significant.  This  shows  that  for  these  LAI  models  for  corn, 
fraction  of  growing  season  based  on  GDD  is  better  than  that 
based  on  HGDD,  and  is  a  far  better  time  scale  than  DAE  or 
fraction  of  growing  season  calculated  using  DAE. 

The  curves  of  maximum  depth  of  depletion  vs.  fraction  of 
growing  season,  based  on  the  1974  through  1984  data  sets,  are 
seen  in  Figs.  9  and  10.  These  curves  were  used  to  estimate  the 
effective  rooting  zone  for  use  in  available  soil  water 
corrections  of  crop  coefficients  (see  Materials  and  Methods 
sect  ion)  . 
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Table  6.  Regression  equations  of  relative  leaf  area  index  for  the 
1981  and  1982  corn  data  vs.  fraction  of  growing  season, 
calculated  by  using  three  methods,  and  days  after  emergence 
(DAE).    R    square   and   root    mean    square    error    are   also    shown. 


Time    base 


Equation 


RMSE 


Fraction 

GDD  t=- 0.1 47+2  1 .6X2-58.5X3+58.3X   -21 . 1X5 


0.7  92  0.174 


Fraction  h 

MGDD  Y=-0.072+13.2X   -23.6 XJ+ 10.6 X* 


0  .779  0. 177 


Fraction 
DAE 


Y=-0  .  163  +  13.8X':-2  4.5X;5+1  1  .OX 


0.779  0.177 


DAE  Y=-0. 282+0. 02 4X-3 . 60 »  10" °X   + 

1.99»10"10X5 


0.635  0.227 
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Table  7.  Regression  equations  of  relative  leaf  area  index  for  the 
1981  and  1982  soybean  data  vs.  fraction  of  growing  season, 
calculated  by  using  two  methods,  and  days  after  emergence  (DAE). 
R    square    and   root   mean   square    error    are    also    shown. 


Time    base 


Eq  uation 


RMSE 


Fraction 
PTU 


Y=0. 126-13. 3X2+63.0X3-82.6X4+32.95  0.901  0.120 


Fraction 
DAE 


Y=-0. 0  48+2  1 . 3X3-38. 7 4+ 1 7 . 6 5 


0.924  0.104 


DAE 


I»0.  1  72-0.0  26X+0.00  1X2-1.  40  •10"7XJt 

+7.28»10"10X5  0.912  0.113 
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Table  8.  Regression  equations  of  predicted  vs.  measured  relative 
leaf  area  index  for  the  1983  corn  data,  predicted  by  using 
curves  based  on  fraction  of  growing  season,  calculated  by  using 
three  methods,  and  day3  after  emergence  (DAE).  R  squares,  root 
mean  square  errors,  and  siginif icance  levels  of  the  F  statistics 
for    the    slope    and   intercept    are    also    shown. 


Time    base 


Fraction 
GDD 


Fraction 
MGDD 


Fraction 
DAE 


Eq  uation 


Y=-0. 045+0 .940X 
Y=0. 025+0. 864X 
Y=0. 13  4+0.7  49X 


RMSE 


Prob>F 
Slope 


Prob>F 
Int. 


0.869        0.138        0.316  0.241 

0.874        0.124         0.015  0  .463 

0.775        0.153        0.001  0.003 


DAE 


Y  =  0. 234+0 .624X 


0.636         0.179         0.0001 


0.0001 
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Table  9.  Regression  equations  of  predicted  vs.  measured  relative 
leaf  area  index  for  the  1983  soybean  data,  predicted  by  using 
curves  based  on  fraction  of  growing  season,  calculated  by  using 
two  methods  ,  and  days  after  emergence  (DAE).  R  square,  root 
mean  square  error,  and  significance  levels  of  the  F  statistics 
for    the    slope    and   intercept    are    also    shown. 


Time    base 


Equation 


RMSE 


Prob>F 
SI  ope 


Fraction 

PTU  Y  =  0. 004+0. 892X 

Fraction 

DAE  Y=0. 028+0. 905X 


DAE 


Y=0. 046+0. 864X 


Prob>F 
Int. 


0.696  0.198  0.219  0.937 

0.784  0.160  0.183  0.509 

0.695  0. 193  0.114  0.364 
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Fig.  9.  Regression  curve  of  maximum  depth  of  depletion 
for  corn  vs.  fraction  of  growing  season  based  on  grow- 
ing degree  days.   The  curve  was  developed  by  using  the 
1974  through  1984  data. 


74 


ai  oo- 

m    : 

A 

x    : 

I  1800- 

M 

U 

M 

1500- 

0  : 

E 

p 

T  1200- 

H 

0 

F   900- 

D 

E 

P   600- 

L 

E 

T 

1  300- 
0 

N 


U  •  ULYSSES 
E  ■  EUDORA 

n  •  mjiR 


I     I    * 


Y-I90.3>2599.3X  -1035. OX. 

R  SQUARE-0.623 

RHSE-234.8 

N-84 


.0    0.1    0.2    0.3    0.4    0.5    0.6    0-7    0.8    0.9    1.0    1.1 
FRACTION  OF  GR0U1NG  SEASON  USING  PTU 


Fig.  10.  Regression  curve  of  maximum  depth  of  depletion 
for  soybean  vs.  fraction  of  growing  season  based  on 
photothermal  units.   The  curve  was  developed  by  using 
the  1974  through  1984  data. 
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The  basal  crop  coefficient  corves  for  corn  are  shown  in 
Figs.  11  through  14.  Using  fraction  of  growing  season,  % 
available  soil  water,  and  Etr  calculated  by  the  method  nsed  to 
derive  the  particular  curve,  we  used  each  of  these  models  to 
predict  the  1983  and  1984  Et  values.  The  predicted  vs.  actual 
Et  regression  lines  are  shown  in  Figs.  15  through  18.  The  TEST 
procedure  was  used  on  these  regression  lines.  The  intercepts 
of  the  lines  resulting  from  using  the  Jensen-Ha i se ,  the 
modified  J ensen-Hai se ,  and  the  Priestley-Taylor-based  curves 
were  not  significantly  different  from  0  (see  Table  10).  The 
TEST  procedure  would  not  test  hypotheses  about  the  predicted 
vs.  actual  Et  regression  line  resulting  from  the  Kcb  curve 
based  on  the  Penman  equation.  Since  its  intercept  is  closer  to 
0  than  those  of  the  other  predicted  vs.  actual  corn  Et  lines, 
it  is  reasonable  to  assume  that  it  also  is  not  siginificantly 
different  from  0.  However,  as  the  F  statistics  listed  in  Table 
10  indicate,  the  slopes  of  these  predicted  vs.  actual  Et  lines 
are  significantly  lower  than  1.  The  prediction  equations 
tended  to  underestimate  Et,  particularly  during  the  middle 
period  of  the  season  when  actual  Et  rates  were  high.  In 
developing  these  crop  coefficient  curves,  measured  Et  rates  of 
various  treatments  within  a  given  y ear-cul t i var  data  set  were 
used.  These  treatments  may  not  have  affected  water  use  of  the 
crop,  and  perhaps  should  have  been  averaged  together  to  avoid 
the  effects  of  unusually  low  measured  Et  values  resulting  from 
variation  in  the  neutron  probe  readings.   A  reevaluation  of  the 
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Fig.     11.    Regression    curve    of    basal    crop    coefficients    for 
corn,    calculated    by    using    the    Penman    reference    Et    equa- 
tion,   vs.     fraction    of    growing    season    based    on    growing 
degree    days.       The    curve    was    developed    by    using    the    1974 
through    1982    data. 
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Fig.  12.  Regression  curve  of  basal  crop  coefficients  for 
corn,  calculated  by  using  the  Jensen-Haise  reference  Et 
equation,  vs.  fraction  of  growing  season  based  on  grow- 
ing degree  days.   The  curve  was  developed  by  usin»  the 
1974  through  1982  data. 
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Fig.     13.     Regression    curve    of    basal    crop    coefficients    for 
corn,    calculated    by    using    the    modified    Jensen-Haise 
reference    Et    equation,    vs.     fraction    of    growing    season 
based    on    growing    degree    days.       The    curve    was    developed 
by    using    the    1974    through    1982    data. 
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Fig.  14.  Regression  curve  of  basal  crop  coefficients  for 
corn,  calculated  by  using  the  Priestley-Taylor  reference 
Et  equation,  vs.  fraction  of  growing  season  based  on 
growing  degree  days.   The  curve  was  developed  by  using 
the  1974  through  1982  data. 
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Fig.  15.  Regression  line  of  predicted  corn  evapotran- 
spiration,  predicted  by  using  the  basal  crop  coeffi- 
cient curve  developed  with  the  Penman  reference  Et 
equation,  vs.  measured  Et  for  the  1983  and  1984  data. 
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Fig.  16.  Regression  line  of  predicted  corn  evapotran- 
spiration,  predicted  by  using  the  basal  crop  coeffi- 
cient curve  developed  with  the  Jensen-Haise  reference 
Et  equation,  vs.  measured  Et  for  the  1983   and  1984  data. 
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Fig.  17.  Regression  line  of  predicted  corn  evapotran- 
spiration,  predicted  by  using  the  basal  crop  coeffi- 
cient curve  developed  with  the  modified  Jensen-Haise 
reference  Et  equation,  vs.  measured  Et  for  the  1983 
and  1984  data. 
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Fig.  18.  Regression  line  of  predicted  corn  evapotran- 
spiration,  predicted  by  using  the  basal  crop  coeffi- 
cient curve  developed  with  the  Priestley-Taylor  refer- 
ence Et  equation,  vs.  measured  Et  for  the  1983 
and  1984  data. 
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Table  10.  Regression  equations  of  predicted  evapotranspiration, 
predicted  by  using  basal  crop  coefficient  curves  developed  with 
four  reference  Et  equations,  vs.  measured  evapotranspiration  for 
the  1983  and  1984  corn  data.  R  squares,  root  mean  square  errors, 
and  the  significance  levels  for  the  F  statistics  for  the  slopes 
and   intercepts    are    also    shown. 


Reference    Et 
eq  uation 


Equation 


RMSE 


Prob>F 
Slope 


Prob>F 
Int. 


Penman 


Y=0 .264+0. 654X         0.737  0.980 


Jensen-Haise      Y=0  .  691 +0  .  688X        0.699  1.133  0.0001  0.065 

Modified 

Jensen-Haise      Y=0 . 53 8+0 . 673X        0.706  1.091  0.0001  0.133 


Priestley- 
Taylor 


Y=0. 520+0. 725X         0.735  1.093  0.0001  0.147 


The      TEST    procedure    in   SAS    did   not    test    hypotheses   about 
this   regression   equation. 
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model  data  and  possible  redevelopment  of  the  curves  is  planned 
in  the  thought  that  a  better  prediction  model  can  be  developed 
from  these  data. 

The  basal  crop  coefficient  curves  for  soybean  are 
presented  in  Figs.  19  through  22,  and  the  regression  lines  of 
the  predicted  Et  rates  using  each  Kcb  curve  vs.  the  actual  Et 
rates  for  the  1983  and  1984  data  are  shown  in  Figs.  23  through 
26.  Table  11  shows  the  F  values  for  these  lines  resulting 
from  the  TEST  procedure.  All  F  values  are  significant  at 
levels  less  than  0.05.  Like  the  corn  curves,  these  soybean  Kcb 
curves  underestimated  Et  rates,  particularly  for  the  higher  Et 
rates  occurring  in  the  middle  portion  of  the  season.  It  is 
thought  that  a  reworking  of  these  data,  in  the  manner  described 
for  corn,  will  improve  the  predictive  ability  of  these  curves. 
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Fig.  19.  Regression  curve  of  basal  crop  coefficients  for 
soybean,  calculated  by  using  the  Penman  reference  Et 
equation,  vs.  fraction  of  growing  season  based  on  photo- 
thermal  units.   The  curve  was  developed  by  using  the 
1974  through  1982  data. 
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Fig.     20.    Regression    curve    of    basal    crop    coefficients    for 
soybean,    calculated    by    using    the    Jensen-Haise    reference 
Et    equation,    vs.     fraction    of    growing    season    based    on 
photo  thermal    units.       The    curve    was    developed    by    using 
the    1974    through    1982    data. 
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Fig.     21.     Regression    curve    of    basal    crop    coefficients    for 
soybean,    calculated    by    using    the    modified    Jensen-Haise 
reference    Et    equation,    vs.     fraction    of    growing    season 
based    on    photothermal    units.       The    curve    was    developed 
by    using    the    1974    through    1982    data. 
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Fig.  22.  Regression  curve  of  basal  crop  coefficients  for 
soybean,  calculated  by  using  the  Priestley-Taylor  refer- 
ence Et  equation,  vs.  fraction  of  growing  season  based 
on  photothermal  units.   The  curve  was  developed  by  using 
the  1974  through  1982  data. 
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Fig.  23.  Regression  line  of  predicted  soybean  evapo- 
transpiration ,  predicted  by  using  the  basal  crop 
coefficient  curve  developed  with  the  Penman  refer- 
ence Et  equation,  vs.  measured  Et  for  the  1983  and 
1984  data. 
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Fig.  24.  Regression  line  of  predicted  soybean  evapo- 
transpiration ,  predicted  by  using  the  basal  crop 
coefficient  curve  developed  with  the  Jensen-Haise 
reference  Et  equation,  vs.  measured  Et  for  the  1983 
and  1984  data. 
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Fig.  25.  Regression  line  of  predicted  soybean  evapo- 
t ranspiracion ,  predicted  by  using  the  basal  crop 
coefficient  curve  developed  with  the  modified  Jensen- 
Haise  reference  Et  equation,  vs.  measured  Et  for  the 
1983  and  1984  data. 
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Fig.  26.  Regression  line  of  predicted  soybean  evapo- 
transpiration,  predicted  by  using  the  basal  crop 
coefficient  curve  developed  with  the  Priestley-Taylor 
reference  Et  equation,  vs.  measured  Et  for  the  1983 
and  1984  data. 
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Table  11.  Regression  equations  of  predicted  evapotr anspiration, 
predicted  by  using  basal  crop  coefficient  curves  developed  with 
four  reference  Et  equations,  vs.  measured  evapotr  anspiration  for 
the  1983  and  1981  soybean  data.  R  squares,  root  mean  square 
errors,  and  the  significance  levels  for  the  F  statistics  for  the 
slopes  and  intercepts  are  also  shown. 


Reference  Et 
equation 


Equation 


RMSE 


Prob>F 
Slope 


Prob>F 
Int. 


Penman  Y=2 .0 81 +0 . 391  X 

Jensen-Haise  Y=  1  .  989+0  .  458X 

Modified 

Jensen-Haise  1=  1 . 9 55+0  .  -i* 44X 


Priestley- 
Taylor 


Y=  1 .  9  1  4+0.425X 


0.278  1.383  0.0001  0.0001 

0.331  1.422  0.0001  0.0001 

0.320  1.424  0.0001  0.0001 

0.264  1.561  0.0001  0.0005 
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Summary    and    Conclusions 

This  study  was  conducted  to  develop  and  test  basal  crop 
coefficient  curves  for  corn  and  soybean  based  on  fraction  of 
growing  season,  calculated  by  using  growing  degree  days  for 
corn  and  photothermal  units  for  soybean.  Four  curves  were 
developed  for  each  crop,  each  one  developed  by  using  a 
different  reference  Et  equation.  The  reference  Et  equations 
used  in  these  curves  were  the  Penman  (Wright  and  Jensen,  1978), 
Jensen-Haise  (1963),  modified  Jensen-Haise  (Burman  et  al., 
1980),  and  Priestley- Taylor  (Kanemasu  et  al.,  1976)  equations. 
These  curves  were  derived  from  water  use  and  meteorological 
data  collected  at  sites  near  Manhattan  and  Tribune,  Kansas 
during  the  years  197  4  through  1982.  The  predictive  ability  of 
the  curves  was  tested  by  using  Manhattan  data  for  the  years 
1983  and  1984.  Other  objectives  were  the  evaluation  of  various 
heat  unit  methods,  the  development  of  curves  of  maximum  depth 
of  depletion  vs.  fraction  of  growing  season  to  estimate 
effective  rooting  depth,  and  the  development  and  testing  of 
relative  leaf  area  index  curves  based  on  fraction  of  growing 
season. 

A      comparison   of    the    coefficients   of   variation   for    fraction 
of    growing   season    to    silking   in   corn,       calculated    by    using   five 
heat    unit   methods    and    days   after    emergence,       led    to    the      choice 
of    the    growing   degree    day    method    (Neild   and    Seeley,       1977).         A 
similar      analysis      of    fraction   of   growing   season      to      beginning 
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bloom  in  soybean  led  to  the  choice  of  maturity  group  II 
coefficients  for  flowering  to  physiological  maturity  for  use  in 
a  photothermal  unit  equation  for  soybean  (Major  et  al.,  1975b). 
Curves  of  corn  relative  leaf  area  index  vs.  four  time 
scales  were  used  to  predict  relative  leaf  area  index  for  the 
1983  data.  F  tests  performed  on  regression  lines  of  predicted 
vs.  measured  relative  LAI  showed  fraction  of  growing  season 
based  on  growing  degree  days  to  be  a  better  time  scale  for  this 
purpose  than  that  based  on  modified  growing  degree  days,  and 
far  better  than  days  after  emergence  and  fraction  of  growing 
season  based  on  DAE.  For  soybean,  curves  of  relative  LAI  vs. 
fraction  of  growing  season  based  on  photothermal  units  and  days 
after  emergence,  as  well  as  one  vs.  days  after  emergence, 
performed   equally   well  in  predicting  the  1983   relative   LAI. 

The  basal  crop  coefficient  curves  for  corn  were  used  to 
predict  the  Et  rates  for  the  1983  and  1984  data.  Linear 
regressions  were  performed  on  the  predicted  vs.  measured  Et 
values.  F  tests  performed  on  these  regression  lines  revealed 
that,  although  the  intercepts  were  not  significantly  different 
from  0,  the  slopes  of  these  lines  were  significantly  lower  than 
1.  This  shows  that  these  basal  crop  coefficient  curves  for 
corn  underestimated  Et.  Likewise,  F  tests  performed  on 
regression  lines  of  Et  predicted  by  using  the  soybean  basal 
crop  coefficient  curves  vs.  measured  Et  for  the  1983  and  1984 
soybean  data  revealed  that  these  curves  underestimated  Et. 

The  basal  crop  coefficient  curves  underestimated  Et  rates, 
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particularly  during  the  middle  of  the  growing  season  when  Et 
rates  were  high.  A  reevaluation  and  possible  redevelopment  of 
the  curves  is  planned.  Averaging  of  Et  rates  of  treatments 
which  did  not  affect  water  use  may  eliminate  some  unusually  low 
outlier  points  in  the  model  data,  thus  allowing  the  curves  to 
reach  higher      maximum      values      at        midseason.  Such        a 

redevelopment  of  the  curves  would  be  particularly  worthwhile 
for  corn,  since  these  curves  were  often  quite  successful  in 
predicting  Et. 
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Table  1  A.  Years,  locations,  soils,  cultivars  (all  Prairie 
Valley  cultivars),  maximum  depths  of  neutron  probe  readings, 
and  references  for  data  sets  used  to  derive  corn  crop 
coefficient    curves. 


Year 

Loc. 

Soil 
(SIL) 

Cultlvar 

Max. 

(  cm) 

Reference 

197  4 

Man. 

Muir 

PV82S 

160 

Stone    et    al .  ,    1  978 

1975 

Man. 

Muir 

PV82S 

160 

Stone    et    al. ,    1 978 

1976 

Man. 

Muir 

PV82S 

160 

Stone    et    al. ,    1 978 

1976 

Man. 

Eudor a 

PV82S 

160 

Anderson    et    al . , 
1982 

1977 

Man. 

Eudora 

PV82S 

160 

Anderson   et    al .  , 
1982 

1978 

Man. 

Eudora 

PV76S 

160 

Anderson   e t    al.  , 
1982 

1979 

Man. 

Muir 

PV76S 

160 

Kufimfutu,     1981 

1980 

Man. 

Muir 

PV76S 

160 

Kufimfutu,     1981 

1980 

Man. 

Eudora 

PV76S 

160 

Kufimfutu,     1981 

1981 

Man. 

Muir 

PV76S 

312 

Hattendorf,     1982 

1981 

Trib. 

Ulysses 

PV76S 

312 

Hattendorf,    1982 

1981 

Trib. 

Ulysses 

PV76S 

312 

Stone,     1982 

1982 

Man. 

Eudora 

PV7  6S 

312 

Redelfs,     1983 

1982 

Trib. 

Oly sses 

PV76S 

312 

Redelfs,     1983 

1982 

Trib. 

Ulysses 

PV7  6S 

312 

Stone,     1982 
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Table  2A.  Years,  locations,  soils,  cultivars,  maturity  groups, 
maximum  depths  of  neutron  probe  readings,  and  references  for 
data    sets    used    to    derive    soybean    crop    coefficient    curves. 


Year 

Loc. 

Soil 
(SiL) 

Cul  tivar 

MG 

Max. 

( cm) 

Reference 

1974 

Han. 

Muir 

Williams 

III 

160 

Mayaki 
1976 

et  al.  , 

1975 

Man. 

Muir 

Amsoy 

II 

160 

Curley , 

1981 

1975 

Man. 

Muir 

Bonus 

IV 

160 

Curl ey , 

1981 

1975 

Man. 

Muir 

Calland 

III 

160 

Curley, 

1981 

1975 

Man. 

Muir 

Clark  Det. 

IV 

160 

Curley, 

1981 

1975 

Man. 

Muir 

Clark  63 

IV 

160 

Curl ey , 

1981 

1975 

Man. 

Muir 

Columbus 

IV 

160 

Curley, 

1981 

1975 

Man. 

Muir 

Pomona 

IV 

160 

Curley, 

1981 

1975 

Man. 

Muir 

Williams 

III 

160 

Curley, 

1981 

1975 

Man. 

Muir 

Woodwor  th 

III 

160 

Curley, 

1981 

1978 

Man. 

Eudor a 

Pomona 

IV 

160 

Stone, 

1978* 

1979 

Man. 

Eudor a 

Pomona 

IV 

160 

Stone  et  al .  , 
1985 

1981 

Man. 

Muir 

Cumberl  and 

III 

312 

Hattend 
1982 

orf , 

1981 

Trib. 

Ulysses 

Cumberl  and 

III 

312 

Hattendorf , 
1982 

19  82 

Man. 

Eudora 

Cumberland 

III 

312 

Redelf s 

,  1983 

1982 

Trib. 

Ulysses 

Cumberl  and 

III 

312 

Redelf s 

,  1983 

Unpublished      work      by      L.  R.       Stone,       Kansas      State      University, 
Manhattan. 
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Table  3A.  Fraction  of  growing  season  to  silking  in  corn 
calculated  by  using  the  growing  degree  day,  modified  growing 
degree  day,  Gilmore  and  Rogers  heat  stress,  Cross  and  Zuber 
heat  stress,  CORNF,  and  days  after  emergence  methods. 


Year 

Cul  tivar 

GDD 

MGDD 

Stress 
GR 

Stress 
CZ 

CORNF 

Days 

1971 

PV82S 

0.521 

0.518 

0.517 

0.517 

0.5  40 

0.518 

1975 

PV82S 

0.456 

0.476 

0.498 

0.512 

0.490 

0.530 

1976 

PV82S 

0.45  7 

0.476 

0.500 

0.510 

0.502 

0  .521 

1977 

PV82S 

0.  472 

0.486 

0.502 

0.511 

0.497 

0.522 

1978 

PV76S 

0.522 

0.529 

0.537 

0.5  42 

0.549 

0.548 

1979 

PV76S 

0.502 

0.511 

0.517 

0.534 

0.530 

0.555 

1980 

PV76S 

0.518 

0.527 

0.5  42 

0.543 

0.547 

0.5  45 

1981M 

PV76S* 

0.477 

0.469 

0.45  9 

0.45  4 

0.494 

0.449 

1982M 

PV76S 

0.460 

0.464 

0.467 

0.480 

0.4  90 

0.495 

1 981  T 

PV76S 

0.572 

0.556 

0.538 

0.534 

0.589 

0.529 

1982T 

PV76S 

0.563 

0.564 

0.563 

0.570 

0.599 

0.576 

1983 

80  5 

0.433 

0.473 

0.516 

0.538 

0.481 

0.565 

1984 

805 

0.464 

0.491 

0.520 

0.527 

0.498 

0.536 

1983 

872 

0.485 

0.526 

0.573 

0.590 

0.539 

0.609 

1981 

874 

0.558 

0.5  80 

0.604 

0.607 

0.595 

0.611 

1983 

980 

0.486 

0.525 

0.570 

0.585 

0.541 

0.602 

1984 

980 

0.551 

0.574 

0.599 

0.601 

0.593 

0.602 

1984 

G4507* 

0.502 

0.525 

0.551 

0.551 

0.548 

0.551 

M  signifies  Manhattan,  T  signifies  Tribune. 


Cultivar  Funks  G-4507.  (These  data  were  not  used  in 
development  of  depletion  curves  and  testing  of  crop 
coefficient  curves.  ) 
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Table  4A.  Fraction  of  growing  season  to  beginning  bloom  in 
soybean  calculated  by  using  the  flowering  to  physiological 
maturity  (F-PM)  coefficients  specific  to  each  maturity  group, 
MG    II    coefficients   for    F-PM,     and    days   after    emergence. 


Year 

Cultivar 

MG 

Specific 
coef  f  s. 
F-PM 

MG    II 
coef  f  s. 
F-PM 

Days 

1974 

Williams 

III 

0.510 

0.446 

0.308 

1975 

Amsoy 

II 

0.554 

0.554 

0.373 

1975 

Bonus 

IV 

0.504 

0.410 

0.310 

1975 

Calland 

III 

0.533 

0.42  8 

0.298 

1975 

Clark  Det. 

IV 

0.502 

0.406 

0.303 

1975 

Clark   63 

IV 

0.50  4 

0.410 

0.310 

1975 

Columbus 

IV 

0.5  90 

0.495 

0.368 

1975 

Pomona 

IV 

0.596 

0.501 

0.377 

1975 

Williams 

III 

0.557 

0.463 

0.324 

1975 

Woodwor  th 

III 

0.625 

0.514 

0.364 

1978 

Pomona 

IV 

0.5  90 

0.457 

0.327 

1979 

Pomona 

IV 

0.512 

0.454 

0.357 

1981M* 

Cumberl  and 

III 

0.503 

0.466 

0.321 

1  981  T 

Cumberland 

III 

0.461 

0.452 

0.310 

19  82M 

Cumberland 

III 

0.531 

0.471 

0.355 

1982T 

Cumberl  and 

III 

0.4  45 

0.417 

0.310 

1983 

Amsoy 

II 

0.509 

0.509 

0.337 

1983 

Union 

IV 

0.674 

0.479 

0.367 

1983 

Bay 

V 

0.460 

0.492 

0.4  92 

1984 

Amsoy 

II 

0.534 

0.534 

0.356 

1984 

Union 

IV 

0  .5  48 

0.463 

0.356 

1984 

Bay 

V 

0.478 

0.46  1 

0.453 

M  signifies  Manhattan,  T  signifies  Tribune. 
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.._„„*.„„  wi  Bi»"iuj  aeaaon  to  tasseiing  and  bl 

stage  in  corn,  calculated  by  using  growing  degree  days  and 
after  emergence.    This  analysis  was  performed  without  the 
and  1984  data.  N=12. 

ister 
days 
1983 

Tasseling                     Blister 

Method    Fraction    Std.  dev.    CV      Fraction   Std  dev. 

CV 

(i) 

ODD  0.450  0.040  8.78 

DAE  0.485  0.032  6.57 


0.6  40  0.0  88  13.66 

0.650  0.078  12.06 
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Table  6A.  Fraction  of  growing  season  to  tasseling  and  blister 
stage  in  corn,  calculated  by  using  growing  degree  days  and 
days  after  emergence.  The  1 983  and  1984  data  are  included  in 
this    analysis.     N=  1  8 . 


Tasseling 


Method      Fraction      Std.    dev.       CV 


Blister 


Fraction      Std.    dev.         CV 


GDD  0.448  0.041  9.08  0.632  0.073 

DAE  0.503  0.042  8.36  0.663  0.066 


11.62 
9.96 


Ill 


Table  7A.  Fraction  of  growing  season  to  beginning  pod  set  and 
beginning  bean  fill  in  soybean,  calculated  by  using 
photothermal  units  and  days  after  emergence.  This  analysis 
was  performed  without  the  7  983  and  1984  data.  N=16  for 
beginning   pod   set.       N=  1  4    for    beginning   bean   fill. 


Beginning 
pod    set 


Beginning 
bean    fill 


Method      Fraction      Std.    dev. 


CV 
(?) 


Fraction      Std.     dev.  CV 


PTU  0.55  9  0.0  42 

DAE  0.481  0.036 


7.51  0.660 

7.44  0.619 


0.049 
0  .0  47 


7.49 
7.65 
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Table  8A.  Fraction  of  growing  season  to  beginning  pod  set  and 
beginning  bean  fill  in  soybean,  calculated  by  using 
photothermal  units  and  days  after  emergence.  This  analysis 
includes  the  1983  and  1984  data.  N=22  for  beginning  pod  set. 
N=20    for    beginning   bean   fill. 


Beginning 
pod    set 


Method      Fraction      Std.     dev. 


Beginning 
bean   fill 


CV  Fraction      Std.     dev.  CV 

(SO  (%) 


PTU  0.568  0.01(1 

DAE  0.497  0.0  48 


7.67  0.671 

9.68  0.636 


0.049 
0.0  49 


7.33 
7.65 
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Table      9A.         Corn      growth      stages    (according   to      Ritchie      and 
Hanway,     1982)    and   dates   of    observance    in    1983. 


Cu 

ltivar 

Stage 

number 

Description 

305 

872 

980 

Planting 

6 

Hay 

6 

May 

6 

May 

VE 

Emergence 

22 

May 

22 

May 

22 

May 

V2 

2    Leaf 

27 

May 

28 

May 

28 

May 

V4 

4   Leaf 

5 

June 

5 

June 

6 

June 

V6 

6    Leaf 

15 

June 

16 

June 

16 

June 

V8 

8    Leaf 

25 

June 

24 

June 

24 

June 

VT 

Tasseling 

7 

July 

15 

July 

18 

July 

81 

Silking 

9 

July 

17 

July 

20 

July 

R2 

Blister 

18 

July 

26 

July 

29 

July 

R3 

Milk 

21 

July 

28 

July 

3 

Aug. 

R4 

Dough 

25 

July 

31 

July 

6 

Aug. 

R5 

Dent 

29 

July 

7 

Aug. 

15 

Aug. 

R6 

Physiological 

maturity 

15 

Aug. 

22 

Aug. 

28 

Aug. 

* 

Tasseling 

(Tassel    visible) 

4 

July 

14 

July 

17 

July 

Growth    stage    in   addition   to    those    of    Ritchie   and   Hanway 
(1982). 
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Table   10A.    Soybean  growth  stages    (according  to   Hanway 
and  Thompson,  1971)  and  dates  of  observance  in  1983. 


Cul  tivar 


Stage 

number 

Description 

Amsoy 

Planting 

16    May 

VE 

Emergence 

26    May 

VO 

Unifoliolate 

leaves 

emerged 

31    May 

V1 

First    trifoll 

ol  ate 

leaves 

13    Jum 

Union 


Bay 


V2 


Trifoliolate  leaves 
at  4  nodes 


21  June 


V3 

Trifoliolate    leaves 

at    6    nodes 

• 

R4 

Beginning    bloom 

26 

June 

R5 

Full    bloom 

8 

July 

R5.5 

Beginning    pod 

development 

13 

July 

R6 

Upper    pod    development 

18 

July 

R7 

Beginning    bean   fill 

24 

July 

R8 

Upper    bean   fill 

31 

July 

R9 

Beans    full    size 

18 

Aug. 

R10 

Physiological 

16  May  16  May 

26  May  26  May 

31  May  31  May 

15  June  14  June 

22  June  22  June 

28  June  30  June 

5  July  28  July 

16  July  4  Aug. 

20  July  15  Aug. 

27  July  18  Aug. 
9  Aug.  23  Aug. 

16  Aug.  27  Aug. 

30  Aug.  10  Sept. 


maturity 


26  Aug. 


12  Sept. 


1  Oct. 


Amsoy  reached  R4  before  this  stage. 


115 


Table    1 1  A.       Corn   growth    stages    (according   to   Hltchie   and   Hanway, 
1982)    and    dates   of    observance    in    1984. 


Cultivar 

Stage 

number 

Descri  ption 

30  5 

874 

980 

Planting 

11 

May 

1  1 

May 

1  1 

May 

VE 

Emergence 

20 

May 

20 

May 

20 

May 

V2 

2    Leaf 
Tasseling 

25 

May 

25 

May 

25 

May 

(Tassel    visible) 

29 

June 

1  1 

July 

16 

July 

R1 

Silking 

4 

July 

17 

July 

21 

July 

R2 

Blister 

17 

July 

23 

July 

30 

July 

R3 

Milk 

21 

July 

29 

July 

4 

Aug. 

R4 

Dough 

25 

July 

3 

Aug. 

10 

Aug. 

R5 

Dent 

29 

July 

9 

Aug. 

18 

Aug. 

F.6 

Physiological 

maturity 

12 

Aug. 

23 

Aug. 

31 

Aug. 
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Table   12A.    Soybean   growth   stages  (according   to   Hanway   and 
Thompson,  1971)  and  dates  of  observance  in  1984. 


- 

Cul  tivar 

Stage 

number 

Descri  ption 

Amsoy 

Union 

Bay 

Planting 

18    May 

18    May 

18    May 

VE 

Emergence 

2  8    May 

2  8    May 

2  8    May 

VO 

Unif  oliolate 

leaves   emerged 

2    June 

2    June 

2    June 

V1       First  trifoliolate 

leaves  12  June      12  June      12  June 


18  June  18  June 

3  July  25  July 
2  Aug. 

20  July  9  Aug. 

28  July  16  Aug. 

4  Aug.  22  Aug. 
14  Aug.  27  Aug. 

maturity                23  Aug.       6  Sept.  3  Oct. 


V2 

Trifoliolate    leaves 

at    4    nodes 

19 

June 

R4 

Beginning   bloom 

28 

June 

R5 

Full    boom 



R5.5 

Beginning   pod 

development 

10 

July 

R6 

Upper    pod    devel 

opment 

20 

July 

R7 

Beginning   bean 

fill 

27 

July 

R8 

Upper    bean   fill 

2 

Aug. 

R10 

Physiological 

Table    13A.       Corn  yield    data   for    1983. 


Ii7 


Cultivar  Plot 


Corn   grain  yield* 


Plant 
po  pul  ation 


805 


872 


9  80 


1 
4 
8 
3 
5 

r 

2 
6 
9 


Kg/ha 
6  548 
6  268 
6  560 
6  452 
6  384 
6  059 
6  021 
4  937 
6    116 


Plants/ha 

69  554 

69  554 
75  459 
73  491 

70  210 
73  491 
73  491 
66  273 
70  866 


Reported   at    15.5it    moisture. 
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Table    14A.       Soybean   yield    data    for    1983. 


Cul  tivar 


Amsoy 


Onion 


Bay 


Plot 


Bean   yield* 


Kg/ha 

1  697 

2  105 
1  691 
1  361 
1  172 
1  80  7 


Plant 
population 


Plants/ha 

355  206 

353  018 

199  913 

255  468 

42  0  385 

378  828 


Reported   at    13.0$    moisture. 


Table    15A.       Corn   yield    data    for    1981. 
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Cultivar  Plot 


Corn   grain   yield* 


Plant 
population 


805 


871 


980 


Kg/ha 

2  964 

3  696 

4  567 

1  822 
4    377 

2  088 
814 

2  596 

3  293 


Plants/ha 

62  336 
57    087 

59  711 
64  691 
77  42  8 
74    147 

60  367 
64   304 

63  648 


Reported   at    15.5?    moisture. 
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Table    16A.       Soybean   yield    data    for    1984. 


Cultivar  pi0t 


Amsoy 


Union 


Bay 


Bean   yield* 


Plant 
population 


Kg/ ha 

Plants/ha 

2 

1  13  4 

189  633 

4 

1  667 

198  163 

9 

1  396 

199  475 

3 

1  266 

146  325 

5 

1  188 

167  323 

7 

1  261 

180  446 

1 

1  217 

121  391 

6 

918 

187  008 

8 

1  450 

162  730 

Reported   at    13.0$    moisture. 
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Table    17A.         Leaf    area    index,       wet   weights,    and    dry    weights   of   805 
corn    oultivar    in    1983.      The    sampling   unit   was   4    plants. 


Date 

LAI 

Wet   wt. 
whole    j 

?1  ant 

Dry    wt. 
leaves 

Dry    wt 
stems 

Dry   wt. 

reproduc- 
tive   parts 

grams  

2    June 

0.05 

8. 

10 

0.77 

0.08 

0 

6    June 

0.04 

17, 

,40 

0.87 

0.44 

0 

14   June 

0.60 

176. 

49 

11  .21 

6.39 

0 

20    June 

1.85 

533. 

,48 

36.84 

22.07 

0 

24   June 

2.78 

1227 

89.99 

62.33 

0.43 

30    June 

2.83 

16  45 

70.5  8 

91.64 

11.90 

8    July 

3.38 

195  9 

91.59 

181.85 

26.87 

13    July 

2.93 

2173 

88.89 

203.71 

48.64 

22    July 

2.77 

2641 

109.45 

206.06 

235.56 

27    July 

2.73 

2875 

108.41 

244.79 

272.92 

4    Aug. 

2.74 

2583 

109.85 

217.18 

289.99 

10    Aug. 

2.18 

2472 

108.95 

228.57 

575.64 

19    Aug. 

0.20 

1716 

118.62 

274.09 

507.62 

29    Aug. 

0 

1039 

54.26 

185  .96 

3  48.31 
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Table    18A.         Leaf    area    Index,       wet   weights,     and    dry    weights   of   872 
corn    cultivar    in   1983.       The    sampling   unit   was   4    plants. 


Date 

LAI 

Wet  wt 
whole 

plant 

Dry    wt. 
leaves 

Dry   wt. 
stems 

Dry    wt. 
reproduc- 
tive   parts 

gr  am  s 

2    June 

0.05 

7 

.50 

0.60 

0.11 

0 

6    June 

0.04 

13 

.50 

1.03 

0.31 

0 

1  4   June 

0.31 

72. 

.07 

5.63 

1.99 

0 

20    June 

1.24 

312, 

.98 

25.39 

12.90 

0 

24   June 

2.02 

673 

49.68 

30.77 

0 

30    June 

3.77 

1644 

87.  12 

67.31 

0 

8    July 

4.94 

2162 

128.31 

141 .68 

8.49 

13    July 

4.24 

2217 

122.41 

141  .81 

12.38 

22    July 

4.59 

306  1 

160.35 

257.49 

133.31 

27    July 

4.49 

3013 

164.00 

313.92 

120  .  1  4 

4    Aug. 

4.21 

3086 

150.86 

292.18 

345.97 

10    Aug. 

3.82 

2370 

144.85 

229.30 

380.67 

19    Aug. 

0.35 

2349 

178.82 

256.46 

673.92 

29    Aug. 

0 

1384 

152 .62 

270.61 

607.37 
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Table    1 9 A.         Leaf    area    index,       wet   weights,     and   dry   weights   of   980 
corn    oultivar    in    1983.    The    sampling   unit   was    4    plants. 


Date 

LAI 

Wet   wt. 
whole    ; 

plant 

D  ry    w  t . 
leaves 

Dry    wt. 
stems 

Dry   wt. 
reproduc- 
tive   parts 

2    June 

0.05 

8. 

.86 

0.69 

gr  am  s 

0.20 

0 

6    June 

0.07 

28. 

,97 

1.99 

0.69 

0 

1 4   June 

0.36 

88 

,45 

6.62 

2.08 

0 

20    June 

1.76 

517, 

.28 

34.71 

18.17 

0 

2  4    June 

2.23 

853 

58.22 

34.54 

0 

30    June 

3.26 

1611 

76.48 

57.88 

0 

8    July 

5.13 

2634 

140.0  3 

144.63 

2.51 

13    July 

5.61 

3752 

159.97 

195  .80 

9.06 

22    July 

5.15 

3974 

191  .28 

388.87 

69.50 

27    July 

3.99 

2  80  7 

152.77 

315.05 

53.77 

4    Aug. 

4.74 

4070 

174.19 

423.95 

215.14 

10    Aug. 

4.48 

4205 

209.04 

484.86 

387.49 

19    Aug. 

1.02 

3300 

177.91 

400.36 

561.65 

29    Aug. 

0.27 

2203 

183.54 

344.51 

635.06 
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Table  20A.  Leaf  area  index,  wet  weights,  and  dry  weights  of 
Amsoy  soybean  cultivar  in  1983.  The  sampling  unit  was  1  meter 
of    row. 


Date 

LAI 

Wet   wt, 

whole    ] 

plant 

Dry   wt. 
leaves 

Dry   wt. 
stems 

Dry   wt. 
reproduc- 
tive   parts 

gr  am  s 

6   June 

0.11 

48 

.  1  4 

5.04 

1.54 

0 

14   June 

0.44 

126. 

.  18 

16.79 

6.75 

0 

20    June 

1  .  12 

250. 

.98 

25.26 

13.60 

0 

22    June 

1.02 

249. 

.28 

25.79 

14.33 

0 

30    June 

2.66 

784 

62. 15 

54.31 

0 

8    July 

2.15 

605 

57.52 

6  5.45 

0 

13    July 

3.57 

1083 

105.62 

11 4.30 

0.34 

22    July 

1.89 

6  46 

70.24 

70.88 

5.15 

27    July 

3.52 

1447 

124.60 

156.50 

53.59 

3    Aug. 

2.82 

142  8 

111.11 

141  .06 

68.83 

9    Aug. 

3.49 

1940 

121  .22 

169.97 

126.93 

16    Aug. 

2.90 

1753 

121.04 

180.82 

181.36 

23    Aug. 

1.98 

16  52 

101 .30 

139.26 

196.57 

31    Aug. 

0.04 

392 

50.25 

95.91 

131.26 
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Table  21A.  Leaf  area  index,  wet  weights,  and  dry  weights  of 
Union  soybean  oultivar  in  1983.  The  sampling  unit  was  1 
meter    of    row. 


Date 

LAI 

Wet   wt, 
whole    ; 

pi  ant 

Dry    wt. 
leaves 

Dry    wt. 

stems 

Dry    wt. 
reproduc- 
tive   part 

6    June 

0.16 

69. 

,60 

7.11 

—  gr  am  s 

1.86 

0 

14   June 

0.37 

104, 

,19 

13.47 

6.59 

0 

20    June 

1.  15 

2  40, 

,48 

24.51 

12.58 

0 

22    June 

0.97 

235. 

48 

24.06 

10.72 

0 

30    June 

2.28 

604 

55.08 

40.5  4 

0 

8    July 

3.09 

82  1 

77.  15 

77.34 

0 

13    July 

2.40 

672 

69.38 

63.32 

0 

22    July 

3.50 

1013 

108.84 

108.95 

0.56 

27    July 

2.21 

775 

80.60 

115.43 

0.46 

3    Aug. 

3.74 

1320 

124.26 

129.00 

6.34 

9    Aug. 

4.47 

1577 

138.99 

181.33 

34.83 

16    Aug. 

4.55 

2051 

158.42 

238.03 

73.53 

23    Aug. 

4.30 

2  7  46 

146.55 

278.95 

157.95 

31    Aug. 

2.77 

1651 

106.99 

199.13 

128.08 

7    Sept. 

0.65 

1354 

92.79 

255.18 

230.76 

13    Sept. 

0.23 

865 

57.29 

190.0  4 

162.49 
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Table   22A.         Leaf    area    index,       wet   weights,    and    dry    weights   of   Bay 
soybean    oultivar    in    1983.      The    sampling   unit   was    1    meter    of    row. 


Dry    wt. 
Wet   wt.  Dry    wt.  Dry   wt.       reproduc- 

Date  LAI  whole    plant  leaves  stems  tive    parts 


6 

June 

0 

.16 

14 

June 

0 

.48 

20 

June 

0 

.63 

22 

June 

0 

.70 

30 

June 

2 

.74 

8 

July 

3 

.68 

13 

July 

3 

.65 

22 

July 

2 

.86 

27 

July 

4 

.27 

3 

Aug. 

5 

.83 

9 

Aug. 

7 

.56 

16 

Aug. 

7. 

,16 

23 

Aug. 

6, 

,41 

31 

Aug. 

5. 

.26 

7 

Sept. 

4. 

.01 

13 

Sept. 

2, 

,82 

23 

Sept. 

1  , 

,68 

3 

Oct. 

0. 

,34 

60. 

12 

126. 

.06 

123. 

88 

160. 

88 

706 

878 

929 

779 

140  2 

1910 

2121 

2710 

2399 

2572 

2199 

16  48 

1637 

1588 

gr  am  s 

5.93  1.99  0 

14.04  5.37  0 

12.97  7.03  0 

16.59  8.44  0 
56.58  48.17  0 

84.60  98.72  0 
114.30  92.97  0 

92.41  9  4.53  0 

129.61  167.56  0 

161.90  237.28  0 

174.41  265.44  0 

231.75  389.55  0.51 

177.15  357.02  6.50 

229.24  384.87  48.41 

169.38  368.60  76.50 

144.11  288.48  53.36 

132.02  283.99  92.92 

90.80  352.33  259.55 
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Table  23A.   Rainfall  record  for  the  Manhattan  site  in  1983. 


Day       May< 


June 


July       Aug. 


Sept.      Oct. 


1 
2 

3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
U 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 

Total 


0 

0 

0.! 

0 

0 

0 

0 

0 

0 

0 
17.0 

6.4 
31.8 

0 

0 

0 

0 
20.8 

1.8 

7.4 
31.8 

1.3 

0 

0 

6.4 

0 

0 

0 

0 

0 

0 

125.2 


0 

0 
10.9 

0 

0 

2.3 

0 

0 

0 

0 

7.6 

0 

0 
39.4 
0 

0 
0 

4.6 

0 

0 

0 

0 

0 

0 

0 

0 

4.1 

6.1 

7.6 

0 


82.6 


0 
0 
0 

1.3 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
19.8 

21  .1 


0 
0 
0 
0 
0 
0 

1.; 

0 

0 

0 

0 

0 

0 

1.8 

0 

0 

0 

0 

0 
0 
10.4 
0 

6.1 
0 
0 
0 
0 
0 
0 
0 

1.1 

67.8 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
36. 
0 
0 

3. 
0 
0 
0 
0 
27. 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


67.8 


May    rainfall    amounts   are   from    the   agronomy    farm    record, 
reldLg.3111   gUaSS   WaS    dissasembled   •"•«•    the   last    neutron 
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Table  24A.       Rainfall    record   for    the    Manhattan    site    in    1984. 


Day 

May" 

June 

July 

Aug. 

Sept . 

0ot.+ 

ID  ID 

1 

0.3 

0 

0 

0 

0 

0 

2 

3.0 

0.3 

0 

0 

20  .3 

2.8 

3 

3.8 

0 

0 

0 

0 

0 

4 

0 

0 

0 

1, 

0 

0 

0 

5 

0 

0 

0 

0 

0 

4.1 

6 

0 

0 

5.6 

0 

0 

8.1 

7 

4.3 

52.1 

0 

0 

0 

0 

8 

0 

1  .0 

0.5 

6. 

4 

0 

__ 

9 

0 

97.8 

0 

0 

0 

__ 

10 

0 

0 

0 

0 

0 

__ 

11 

0 

0 

0 

0 

0 

__ 

12 

0 

0 

0 

0 

0 

.. 

13 

0 

14.2 

0 

0 

0 

__ 

14 

0 

0 

0 

0 

109.5 

.. 

15 

6.1 

41  .9 

0 

0 

1.5 

_. 

16 

0 

0 

0 

0 

0 

__ 

17 

0 

2.0 

30.0 

0 

0 

__ 

18 

0 

0 

0 

0 

2.0 

__ 

19 

50  .8 

0 

0 

0 

0 

__ 

20 

0 

3.0 

0 

0 

0 

__ 

21 

0 

24.9 

0 

25. 

4 

0 

__ 

22 

11.7 

20.8 

0 

0 

0 



23 

0 

16.8 

0 

0 

0 

__ 

24 

0 

0 

0 

0 

0 

__ 

25 

2.3 

0 

0 

0 

0 

__ 

26 

15.5 

3.3 

0 

0 

0.8 

__ 

27 

25.4 

2.5 

0 

0 

10.4 

__ 

28 

0.5 

0 

0 

0 

0 

__ 

29 

0 

0 

0 

0 

0 

__ 

30 

0 

0 

0 

0 

0 

__ 

31 

0 

0 

-- 

0 

8 

-- 

-- 

Total 

123.7 

280.6 

36.1 

32. 

144.5 

May    rainfall    amounts   are    from    the   agronomy    farm    records. 


The    rain    guage    was    disassembled    after    the    last       neutron 
probe    reading. 
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i«  'iF";F;<o01r7Hi:,s45FYp^^-3s,,"-d3i76*,"vG--"-= 

'03  AC.15FP*M5FP; 
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TusiL2fA-       SAS    Pr°8ram    f0r    calculating    "feren 
using    four    equations.  n 


ce     Et    by 


1  DATA   STR: 

*  TMINC»(TMINF-32I».J56; 

5  TCa-(TFa-32)».S56i 

»  «RKM«w»MI«l.61; 

KR2-t«KW»(3.28»*.2); 
5  uSf"1    HANHATT4N    ANE«a«6TER    HI*    RUN    TO    A   2     ,ETER    HEIGHT; 

*C  PROC  PRINT; 

»CTE=    THE   PRCCEOURE   prINT   USEO   0.2,    SECONDS    ANO   202K    ANO    PR!NTEO    PAGE    U 

♦1  DATA    TkO: 

*J  SET   ETRj 

*3  CAVG»<TNAXC»TNINC)/2; 

**  L-3?5-0.Sl«CAVG; 

«      s-sslSSiaiSi,""  tri8une  site  «  «»••« 

*8  *«<a.3  86»p|/i.; 

52  ESAve-IESMIN-ESNAX)/2;  ♦O.a072)«»8-0.0000 »•( l.a»TMAXC-*8H-0.0Ol3l6l i 

RH?STa"H:;[oo?0n8*TCa^-a°721**9"-0000"*<^''^«^^0.001316.= 

55  EDP8-ESa»RHFRACT8: 

"  »    N£T    RAOIATICN    EOUATICNS! 

"  *NCO«N1"0-861»S01.RAO- 103.92; 

5|  "NCCRN2»0.84a«SaLBA0-lV».*s; 

"  *NWRN3-0.7t6»S0LRA0-?9.89: 

}«  SMSrfRG-O.725«S0LRAO-5a.7«.; 

*1  RNSOY82»0.3C5«SOLRAO-137.5<>; 

M  RNSORG2-J.a37»SOLRAO-132.34; 

"  EpLT?,f^,Ur.*!ES'JSIN<:   THE    PENMAN    SCUATICN; 

EPENCl.(0ELTA/IOtCTA.r)).RNCORN1.(r/IoEUTA;y,1.15.36.(o.75.o.OU5.-R2).(ESAVE-Co 

J.  P8Pff2-,D£^^'O^^"*^CU»N2.(V/,0EtTA.rnn5.36.(0.75.o.=  uS.WR2I.(ESAVE-^ 

^  6PENC3.,0ELTA/(CELTA>V,,.RNCORN3»,V/(UEl.TA.nJ.l5.36.(0.75.O.OU5.WR21.lESAVE-E0 

=  p:fr-,CELrA/'QELr4-',«N^^<-'°^-n.I5.36.,o.7„o.ou5-,R21.„SAvE-Eo 

2  P9?f2-1DEL^'-^-'-NSOVa2Mr/,UELTA.vn.15.36.<O.75.0.0U5.WR2).(ESAvE-E0 
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2  S    A   S        L    0  G           OS   SAS    32.3                       OS/HVT   JOB    XPRS93!7    STEP    SAS                PROC 

IS  ^f*G2""5ELTA/,DELT4*1"",RNS0R«<-|Y/"'ELTA»rl,.15.34«(0.75»0.aU5.v<R21.,ESAV6-ED 

74  LflN-L/10: 

t!  =.!I5,CAU:iJLATED    USING    me    PRIESTLEY- TAYLOR    EQUATION: 

78  SPRrYCl-U32«(DELTA/(0ELTA*r)j«RNC0RNi; 

"  EPRTYC2-L32«<DELTA/tOELTA*Y)  l«RNCnRN2; 

80  6|,RTYC3-l.32-(DELTA/(0ELTA*YJ)«RNCORN3; 

81  EPTS0YR.1.32*(0ELTA/(CELTA*Y))«RNSOYRG; 

82  EPTS0Y2-L32*<0ELTA/(0ELTA<-YI)»RNSaYS2: 
«  E  PTSRG2"! .32* { DELTA/ ( OE  LTA+Y ) ) *RNS0RG2 : 
8+  FAVG-<TP*XF»TMINF>/2; 

11  *    ETR   CALCULATED*  USING    THE    JENSEN-HAISE    EOUATION! 

8t  JHIY»(0.014»FAVG-0.37)»SDLRA0: 

8T  TCWNAX»33.2s 

88  TCli)1IN-20.0; 

a  «sr»j^^^T^a^,OTBis5,'"  july  at  ^"^  ™ese 

95  CH-50/1EW1AX-EMMINJ; 

94  CT»1/(C1+7.3«CH); 

97  TX— 2.5-0.  l*»<EKMAX-€WHIN)-(ELEVM/550); 

'8  «0OJHLY=<T»(CAVG-TXI«S0LRAO! 

TOO  IttSSKES?   ETR    mU"   FR0"   UNGUYS    T°   ""    PE"    0AY*      PREfI*    «   SIGNIFIES 

IC1  PENCCRl»ePEM:i./LWK: 

102  IF    P6NC0RK0    THEN    PENCORl-0; 

1C3  ACPENCI+PENC0R1; 

104  PENCCR2-EPENC2/IMH; 

105  IF    PENC0R2<0    THEN    PENC0R2-0: 
104  ACPENC2»PENC0R2: 

107  P6NCCR3-EPENC3/LMM; 

lcs  IF    PENCCR3<0    THEN    PENC0R3-0; 

109  ACPENC3-»PENCGK3; 

no  pensoyr.epsoyr/lhn: 

hi  if  pensoyr<o  then  pensoyr-o: 

112  acpensry*pensoyr; 

113  penscy2-epscy2/lnn: 

u*  if  p£nsoy2<0  then  pensoy2»0: 

us  acpensy2*penscy2: 

u4  pensrg2-epsrg2/lhm; 

Jl'  IF    PENSRG2<0   THEN    PENSRG2.0: 

118  ACPENSG2*PENSRG2; 

U9  PTCORNI-EPRTYCl/LMN: 

<•*?  IF    PTCORNKO    THEN    PTCCRM-O: 

121  ACPTCORUPTCQRNl: 

122  PTCGRN2-EPRTYC2/LJW 

123  IF    PTC0RN2<0    THEN    PTCCRN2-0; 
1-2+  ACPTC0R2»PTOJRN2: 

125  PTCCRN3»SPRTYC3/L«C: 

124  if    PTCORN3<0    THEN    PTC0RN3-0-. 

127  ACPTC0R3+PTC0RN3: 

128  PTSCYRG-tPTSOYR/LHI": 

129  ip    PTSOYRG<0   THEN    PTSCYRG-O: 
120  ACPTSYRGtPTSOYRG; 

131  l>TSOYB2»E?TSOY2/Lmi; 
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3  SAS        L    0  i  OS    SAS    82.3  OS/MVT    JOB    XPRS93I7    STEP    SAS  PROC 

132  IF   PTS0YB2<0    THEN    PTSOYB2=0; 

133  ACPTSYB2«-PTS0Y82: 
13*  PTSRGM2-EPTSRG2/LHM; 

135  IF    PTSRGH2<0    THEN    PTSRGM2»0; 

136  ACPTSRG2*PTSRGK2: 

137  JHHB-JHLY/LMM; 

138  ACJHKM+JHHMi 

139  MOOJHMM=HODJHLY/LMN; 
1*0  ACHCDJH+MOOJH/1M; 

NCTE:  DATA    SET    MURK. TWO    HAS    30   OBSERVATIONS    AND    83    VARIABLES.    28    OBS/TRK. 

NCTE:  THE    DATA    STATEMENT    USEO   0.81    SECCNOS    ANO    208K. 

1*1  PROC    PRINT; 

1*2  VAR    OATE  THAXF   TM1NF   RH3   TC8    KRMI    TMAXC    THINC   WRKH    WR2    SOLRAO 

1*3  DATE  CAVG    L    ELEVC   P    Y   DELTA    ESCIN    ESHAX    ESAVE    ES8   RHFRACT8    EDP8 

1**  DATE   RNC0RN1    RNC0RN2    RNC0RN3    RNSOYRG   RNS0YB2   RNS0RG2    EPENC1    EPENC2    EPENC3 

1*5  OATE   EPSOYR   EPS0Y2    EPSRG2    LMM    EPRTYC1    EPRTYC2    EPRTYC3    EPTSOYR    EPTSUY2    EPTSRG2 

1*6  DATE    FAVG   JHLY    TCWMAX    TCWM1N    EWMAX    EWM  IN    CI    CH    CT    TX    MOOJHLY   PENC0R1    ACPEMC1    DAT 

1*7  E  PENC0R2    ACPENC2    PENCDR3    ACPENC3    PENSOYR    ACPENSRY    PENS0Y2    ACPENSY2    CATE 

1*8  PENSRG2    ACP1NSG2    PTC0PN1    ACPTC0R1    PTCURN2    ACPTCCR2    PTCQRN3    ACPTC0R3    PTSOYRG 

1*9  DATE   ACPTSYRG    PTSCYB2    ACPTSYB2    PTSRGH2    ACPTSRG2    JHPM    ACJHMM   MODJHMM    AOTODJH    DATE 

150 

NCTE:  THE    PROCEDURE    PRINT    USED    1.50    SECONDS    AND   220K    ANC    PRINTED    PAGES    2    TO    5. 

NCTE:  SAS    USED    220K   MEMORY. 

NCTE:  SAS    INSTITUTE    IMC. 
SAS    CIRCLE 
PO    SOX    3000 
CARY,    N.C.       27511-8000 
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'strati™    tlil    ll°STam    f°r    eStiraatinS    "ybean    evapotran- 


l  OATA    EST; 

2 


»  fJ&sJiOm""   "FGS    ""    "    PEN"    JH4    ™M    "»■ 


KA- 11.00  USW»UJ/«. 615; 
5  PEN.PEKK/100; 

4  JH-JHA/100; 

N0DJH-M0DJHA/100; 

PT-PTY/100; 
*  ETMn-ET/LOO: 

10  CARDS: 


ss  S'js  SB»^».a,aga,^Iji!"««-  -  ««*. 


PRQC  PRINT; 


"CTE:    THE   PROGECUAE    M»T    USEC   0.,3    SECQf|US    ^    ^    ^    ^^    ^    ^ 

2*  DATA   TWO; 
*'  SET    EST; 

^5"l?5s;prt^^t'F\1I"3*,FGS*faS*FGS'-42-89"27«1*l"S.FGS.PGS.fGS)^C.6S, 

»  PENKCW-PENKC3»KA; 

J'  fTPENW-PENKCH*PEN; 

|0  PENKSl-(O.<J-PS(KC8)»0.a; 

31  ?EJIKS2--(0.9-PEMKCe!*0.5; 

32  PE.,lXS3«(<3.9-PENKCaj»0.3  ; 
32  PENKCM1.PEN^31«PEN; 

JJ  PENK,«M2.PEf,KS2«PEN: 

3-  PENKNH3-FENKS3«PEN: 

"  PEWOr.PENKNMltPEfiKHm.PENKMMa; 

3'  JHKCH-JHKC3»KA; 

40  6Tjm*jHKCW«JH; 

*i  J"KS1-[G.9-jhkC3J»0.3; 

*|  -IHKS2-<C.<)-JHKC3)«0.5; 

*-  jmks3.( a.*-jHKce;*o.3 ; 

**  JHKMNl«jHKjI»*; 

*!  jhkk«2-jhks<»jh; 

♦«  Jhkivu-jhksJ'jh: 

*'  JHKIar.JHK«l1l*JHK,MM2*jHHMH3; 

!C  l'OOKCW-nDDKC3»KA; 

J!  =W00M-«00KC^»H00JH: 

'«  "OCKSl«(0.9-HC;.<C3l-0.3; 

'-  -1CJKS2-  (C.9-1CCKC3)«0  .5  : 

«  fOUK33- (0.9-nODKCB  I  «0  .3  ; 

1  -  "QOKcrti.fnoKs  i»  »\oojh  ; 
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2.  SAS        L    a   G  OS    SAS    82.3  CS/HVT   JOB    XPRS0*5*    STEP    SAS  PRCC 

54  N0DKMM2-HO0KS2*MOCJH; 
57                      MClOKH.13»HaOKS3*MGC.JH; 

5!  .1COXTOT«»«OOKMI11»*C)CKM»2»MaDKKI43: 

55  PTKCS-0. 277 2*393*23. 772 11660* (FGS'FGS'FGSI-vi.dCOOO 63** (FGS*FGS*FGS*FGS 1  +  2 1. 13*9 
6C  8830*<FGS*FGS»FGS»FGS»FGSJ  ; 

61  PTKCX-PTKCa*KAi 

62  ETPTW"PTKCW*PT; 

63  PTKSl-(0.9-PTKC3)»0.a  : 
6*  PTKS2-(C.9-PTKC3)*0.5 ; 

65  PTKS3*(C.9-PTKC3)*C3: 

66  PTKWU»PTKSl*PT: 

67  PTK!W2*PTKS2»PT: 
6!  PTKMK3"PTKS3*PT: 

6?  PTT0T-PTK«H1*PTKMM2+PTKMM3: 

iCTE:    DATA   SET    *0SK.TUC    HAS    12   JESERVATTCNS    AND    56    VARIAflLES-    *2    OBS/TRK. 
NCTE:    THE    OATA    STATEMENT    USES   0.*2    SECCNCS    AIVC    203K. 

7C  PSCC    PRINT: 

NCTE:    THE    PROCECURE   PRINT    uSEO    0.57    SECONOS    AND    222K    ANO    PRINTED    PAGES    2    TO   3. 
NCTE:    SAS    LSEO    222K   rtE,»URY. 

NCTE:    SAS    INSTITUTE    I.'IC. 
SAS    CIRCLE 
PC   3CX    3CC0 
CARY,    N.C.      27511-SOOO 
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ABSTRACT 

Basal  crop  coefficient  (Kcb)  curves  for  corn  (l^a  aay_s  L.  ) 
and  soybean  (Glycine  max  (L.)  Merr. ) ,  based  on  the  time  scale 
fraction  of  growing  season  calculated  by  using  growing  degree 
days  for  corn  and  phot  othermal  units  for  soybean,  were 
developed  from  neutron  probe  water  use  data  collected  at  a  site 
near  Manhattan,  Kansas  during  the  years  1971  through  1982  and 
at  a  site  near  Tribune,  Kansas  during  the  years  1981  and  1982. 
The  crop  coefficients  were  basal  because  they  were  developed 
from  data  collected  during  times  of  minimal  evaporation  and 
were  corrected  for  the  lowering  effects  of  inadequate  soil 
water.  The  time  scale  was  chosen  with  the  thought  that  such 
curves  would  not  be  subject  to  shifts  due  to  cultivar  and 
climate  related  differences  in  crop  development.  An  examination 
of  the  coefficients  of  variation  (CV)  of  fraction  of  growing 
season  to  silking  in  corn  and  beginning  bloom  in  soybean, 
calculated  by  various  heat  unit  methods  as  well  as  days  after 
emergence  (DAE),  led  to  the  choice  of  the  methods  mentioned 
above.  Four  reference  Et  (Etr)  equations  were  used  in  these 
curves.  These  were  the  Penman,  Jensen- Haise ,  modified  Jensen- 
Haise,  and  Priestl  ey- Taylor   equations. 

With  the  inputs  of  fraction  of  growing  season,  Etr  values 
calculated    from    meteorological    data,     and     correction 
coefficients   calculated  from  available  soil  water   percentages 
and   rainfall  amounts,   the  basal  crop  coefficient  curves   were 


used  to  predict  Et  rates  for  data  collected  at  the  Manhattan 
site  in  1983  and  1984.  The  predicted  Et  rates  were  regressed 
on  the  measured  Et  rates  and  F  tests  were  performed  on  the 
linear  equations.  For  corn,  the  intercepts  of  these  lines 
were  not  different  from  0,  but  the  slopes  were  significantly 
lower  than  1.  In  the  case  of  soybean,  the  intercepts  of  these 
lines  were  significantly  greater  than  0  and  the  slopes  were 
significantly  lower  than  1.  This  showed  that  the  basal  crop 
coefficient  curves  underestimated  the  corn  and  soybean  Et  rates 
for  the  1983  and  1984  growing  seasons. 

Curves  of  maximum  depth  of  depletion  vs.  fraction  of 
growing  season  for  corn  and  soybean  were  developed  from  1974 
through  1984  water  use  data  collected  at  the  Manhattan  site  and 
from  1981  and  1982  data  collected  at  the  Tribune  site.  Curves 
of  corn  and  soybean  relative  leaf  area  index  (LAI)  vs.  fraction 
of  growing  season,  calculated  by  various  heat  unit  methods  and 
DAE,  were  developed  from  1981  and  1982  data  collected  at  the 
Manhattan  and  Tribune  sites.  The  equation  based  on  fraction 
of  growing  season,  calculated  by  using  the  GDD  method, 
successfully  predicted  relative  LAI  of  the  1983  data  collected 
at  the  Manhattan  site  while  the  other  methods  failed.  The 
soybean  relative  LAI  models  predicted  the  Manhattan  1983 
relative  LAI  equally  well. 


